Historic, Archive Document
Do not assume content reflects current
scientific knowledge, policies, or practices.
U.S.D.A. Forest Service
Research Paper RM-42
1968
simulating the management of
even-aged timber stands
Rocky Mountain Forest and Range Experiment Station
U.S. Department of Agriculture
Forest Service
PBBHE = Q.9546£*ALaG10<BBHa> - 0. 1Q640*ALDG1G<PRET> + 0,26959
llll
III
II
I I I
I II II
oooooooooooooo|oooooooooooo||oooooooo|oo|ooloooooo||ooo|oooo|ooooooooooooooooooo
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 4$ 47 43 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80
1 1 11 1 1 1 1 11 1 1 1 1 1 1 11 11 1 1 | 11 1 | 1 1 1 1 1 1 11 1 1 1 1 | 1 11 1 1 | 1 1 1 | 1 1 11 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 11 1 1 1 1 11
22222222|22222222222|222222222|2222222222222222222222222222222|22222222222222222
333333333333|33|3333333|33333333333333|3333333|33333333|33333|333333333333333333
4444444|4444444444|44|444444||444|44444444|4|444444|4444|44444444444444444444444
5555555555|555555|555555555555555555555555555555555555|5555555555|55555555555555
66 6 6666666666666666|6666|6666666|66666666|66666|666666666666666|66666&6666666666
777777|777777777777777777|7777777777777777777777|777|777777777777777777777777777
888888888|88|88|88888|888888|88|S|8888|88888|888888|8888|8888|8888888888888M888
9999999999999999|999999999999999999999999999999999999|9999999999|9|9999999999999
' 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 3? 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80
GIOBF NO 1 STANDARD FQrM 5081 PRINTED IN U S A
U.S.D.A. Forest Service
Research Paper RM-42
1968
Simulating the Management of
Even-Aged Timber Stands
By
Clifford A. Myers, Principal Mensurationi st
Rocky Mountain Forest and Range Experiment Station
Central headquarters maintained at Fort Collins in cooperation
with Colorado State University.
CONTENTS
Page
Introduction 1
Description of program MANAGD 2
Data deck 8
Description of test problem 12
Literature cited 15
Appendices:
1. Program MANAGD 16
2. Output of test problem 23
Simulating the Management of Even-Aged Timber Stands
Clifford A, Myers
Introduction
Good decisions that lead to efficient, low-
cost operations are as necessary for timber
management as for any other business activity.
Gathering of information on which to base such
decisions is unusually difficult because of the
peculiar characteristics of the business of grow-
ing wood. The time needed to mature a crop
is long, and, with many tree species, results of
changes appear slowly. Decisions must be made
for large areas of land and for numerous pos-
sible combinations of biological and economic
conditions. A procedure whereby a manager
can examine quickly the probably outcomes of
many variations in management should be of
value. Such a procedure appears especially
useful when the computation time required is
about 10 seconds for each combination of pos-
sible alternatives examined. This approach to
managerial decisionmaking is made possible
by electronic computers and the technique of
simulation.
Simulation involves the creation and opera-
tion of a model that is similar in relationships
to the natural system studied. The system may
be a factory, a forest, or any one of countless
other subjects of interest. Mathematical simu-
lation, where the model is a set of equations,
is a powerful tool of management and investiga-
tion. Mathematical models are used to discover
new facts and to test alternatives. Time can
be compressed so a sequence of events can be
studied in a fraction of the time required for
operation of the natural system, thereby permit-
ting projections into the future (Choraf as 1965).
Manipulation of a model by simulation is
especially useful where: (1) the actual system
should not be changed until decisions have
been made; (2) possible changes are numerous
and are costly in money, effort, and time; and
(3) exact data cannot be obtained because of
cost or complexity of the relationships involved
(Churchman et al. 1957, Dane 1966).
Several authors have presented the prin-
ciples of simulation to foresters, and have shown
how forestry simulators may be constructed.
These authors have discussed the development
of models (Chappelle 1966) and the systems
approach to forestry problems (O'Regan et al.
1966). Others have described simulators of
timber production activities, and presented in-
formation obtained from their operation (Clut-
ter and Bamping 1966, Gould and O'Regan
1965). Thompson's (1966) discussion of the
need for tests of alternatives in forest regulation
points out a value of simulation for research.
Validity of many long-accepted concepts of
timber management can be tested quickly,
cheaply, and for a wide variety of conditions.
Gould and O'Regan (1965) performed a
special service to timber managers by publishing
a simulator at a time when examples were
relatively inaccessible to most foresters. Their
computer program will long serve as a source
of ideas and procedures; it has so served for
the program presented here.
Program MANAGD (appendix 1) is a means
for simulation of the management of even-aged
stands for roundwood and saw logs. It contains
provisions for stand growth, thinnings, harvest
cuts, planting of nonstocked areas, and other
changes in forest conditions. Inputs to the
program permit wide choice in the management
alternatives and stand conditions to be ex-
amined. Possible options and alternatives are
presented in the appropriate sections of the
program description.
The program was written in A.S.A. Fortran
IV (X3.9-1966) and tested on a CDC 6400 compu-
ter. 2 Program organization will permit ready
modification to fit local tree species and com-
puting equipment. Once the program is in use,
improved or additional data and procedures can
be inserted easily, as they become available.
1 Trade names and company names are used
for the benefit of the reader 3and do not imply
endorsement or preferential treatment by the
U. S. Department of Agriculture.
- 1 -
Description of Program MANAGD
Program MANAGD consists of a main pro-
gram and 11 subroutine subprograms (appendix
1). Content and purpose of each routine are
given in the sections that follow. Definitions
of variable names are listed with the source
program in appendix 1 and in the listing of
contents of the data deck. The test problem
described on page 14 and reported in appendix
2 provides additional explanation of the program.
The terms batch, test, and game are used to
identify individual simulations performed with
various groupings of alternatives. The BATCH
name identifies the entire group of tests and
games to be completed as a single job by a
computer. A test consists of one or more
games, all of which are based on a single yield
table and one set of stumpage prices. Games
of a test differ from one another in distribution
of acres by age classes, area planted, and limi-
tations on the annual cut.
Main Program
The main program calls 10 subroutines in
proper sequence, and uses counter UK to call
an eleventh subroutine (OUTPUT1) at specified
intervals. The first four subroutines called
read the data deck, compute and print a yield
table, and compute and print potential volumes
per acre at each year of stand age. A fifth
subroutine generates a working circle with the
specified number of acres in each age class.
Subroutines YEARS and ANUAL are the
dynamic components of the program. They
provide for stand growth, thinnings, harvests,
losses, and other changes in volume and value.
Four subroutines produce a record of results.
Numbers of acres in each age class are printed
at the end of the first year of a game and at
the end of each decade. Volume and dollar
values are printed at the end of each game.
An optional subroutine (SUMRY) prints selected
values from several games of a test together
to simplify comparisons.
The main program enters BATCH name and
the number of tests in a batch into computer
memory.
Subroutine INPUT1
INPUT1 is called once each test to read
values that apply to all games of a test. Values
entered include stumpage prices, minimum
commercial volumes, and items used to compute
a yield table. Controls on the program are
entered as number of games in a test, number
of years in a game, and number of columns
of OUTPUT2 to be summarized by SUMRY.
Values read by subroutine INPUT1 are
printed as part of the output of other sub-
routines.
Subroutine YIELD
Subroutine YIELD reads four tables and
computes and prints yield tables for managed,
even-aged stands. It is called once each test
to produce the yield table that will apply to
all games of the test. Values in each yield
table reflect prior decisions on the frequency
and intensity (table 1) of intermediate cuttings.
Computations performed by the subroutine
follow procedures described in detail elsewhere
(Myers 1966, Myers and Godsey 1968). Average
stand diameter (d.b.h.) and number of trees
per acre just before initial thinning are esti-
mated from measurements of stands of suitable
ages and densities. Number of trees and d.b.h.
are used to compute initial and subsequent
basal areas per acre. Total cubic volumes per
acre are computed from basal area, d.b.h., and
tree height (table 2) by means of a stand
volume equation (statements for TOTO and
TOTT). Total cubic volumes are multiplied
by factors to convert them to merchantable
cubic volumes (table 3) and board feet (table
4).
Increase in d.b.h. due to initial thinning is
computed by the iterative process of FORTRAN
statements 6 to 12, inclusive. The equation
for PDBHE provides an estimate of diameter
after thinning when diameter before thinning
and the percentage of trees retained are known
(Myers 1966). Successive percentages are tested
until d.b.h. after thinning, number of trees,
and basal area agree with the desired thinning
intensity entered as THIN by INPUT1 and as
defined in table 1. Equations for DBHP, such
as statement 7, express that part of table 1
where basal area is less than 66.3 square feet.
Table 1 gives residual basal area after thin-
ning for various average stand diameters. The
values represent one possible series of densities
that could be used to guide successive thinnings
in a stand. The series is labeled growing stock
level 80 to indicate that reserve basal area is
80 square feet per acre when d.b.h. after cutting
- 2 -
Table 1. — Basal areas after intermediate cutting in relation to average stand diameter
Growing stock level 80 for Black Hills ponderosa pine
Average
Basal
Average
Basal
Average
Basal
Average
Basal
stand d.b.h.
area
stand d.b.h.
area
stand d.b.h.
area
stand d.b.h.
area
after cutting
per
after cutting
per
after cutting
per
after cutting
per
( Inches )
acre
( Inches )
acre
(Inches)
acre
(Inches)
acre
Sq . f t .
Sq . f t .
Sq. ft.
Sq. ft.
2.0
12.3
4.0
35.2
6.0
56.6
8.0
72.5
2.1
13.3
4.1
36.4
6.1
57.6
8.1
73.2
2.2
14.4
4.2
37.6
6.2
58.5
8.2
73.8
2.3
15.5
4.3
38.7
6.3
59.4
8.3
74.4
2.4
16.6
4.4
39.8
6.4
60.3
8.4
74.9
2.5
17.7
4.5
41.0
6.5
61.1
8.5
75.4
2.6
18.9
4.6
42.2
6.6
62.0
8.6
75.8
2.7
20.1
4.7
43.4
6.7
62.9
8.7
76.3
2.8
21.3
4.8
44.6
6.8
63.8
8.8
76.7
2.9
22.5
4.9
45.7
6.9
64.7
8.9
77.1
3.0
23.6
5.0
46.7
7.0
65.5
9.0
77.5
3.1
24.8
5.1
47.7
7.1
66.2
9.1
77.9
3.2
26.0
5.2
48.8
7.2
67.0
9.2
78.2
3.3
27.2
5.3
49.9
7.3
67.7
9.3
78.5
3.4
28.4
5.4
50.9
7.4
68.4
9.4
78.8
3.5
29.6
5.5
51.8
7.5
69.1
9.5
79.0
3.6
30.7
5.6
52.8
7.6
69.8
9.6
79.2
3.7
31.8
5.7
53.7
7.7
70.5
9.7
79.5
3.8
32.9
5.8
54.7
7.8
71.2
9.8
79.7
3.9
34.1
5.9
55.6
7.9
71.9
9.9
79.9
10.0+ 80.0
Table 2. — Average height of dominant and codominant trees at
various ages, Black Hills ponderosa pine
Main
Site index
class
stand age
(Years)
40
1 50
| 60
70
- - - - - Feet
10
4.5
4.5
4.5
4.5
20
9
10
12
16
30
11
16
20
26
40
17
22
28
35
50
21
28
35
43
60
26
33
41
50
70
30
38
47
56
80
34
43
52
61
90
37
47
57 ~
66
100
40
50
60
70
110
43
53
63
74
120
45
56
66
77
130
46
59
69
80
140
48
61
71
83
150
50
63
73
86
160
51
64
75
88
170
52
65
77
90
180
53
66
78
91
- 3 -
Table 3. — Factors for conversion of stand volumes in total cubic feet to
merchantable cubic feet per acre, Black Hills ponderosa pine
Average
!
Average
Average
stand
Ratio
of |
stand
Ratio
of !
stand
Ratio of
diameter
merchantable 1
diameter
merchantable j
diameter
merchantab le
(Inches)
to total
volume :
(Inches)
to total
volume |
(Inches)
to total volume
5.0
0.332
8.1
0.849
11.9
0.940
5.1
.355
8.2
.856
12. 1
.941
5.2
. j / /
8.3
. 862
12.4
.942
5.3
.400
8.4
.868
12.7
. 943
5.4
.422
8.5
.872
12.9
.944
5.5
.444
8.6
.876
13. 1
.945
5.6
.465
8.7
.880
13.3
.946
5.7
.487
8.8
.884
13.5
.947
5.8
.508
8.9
.888
13. 7
.948
5.9
.530
9.0
.892
13.9
.949
6.0
.552
9.1
.896
14.2
.950
6 . 1
.575
9.2
.899
14.4
.951
6.2
.597
9.3
. 902
14. 7
.952
6.3
.618
9.4
.906
14.9
.953
6.4
.639
9.5
.910
15.2
.954
6.5
.659
9.6
.913
15.4
.955
6.6
.678
9.7
.916
15.8
.956
6 . 7
. 694
9.8
.920
16.3
.957
6.8
.710
9.9
.923
16.8
. 958
6.9
.725
10.0
.926
17.3
.959
7.0
.740
10.1
.928
17.8
.960
7. 1
.753
10.2
.930
18.3
.961
7.2
.766
10. 3
.931
18.8
.962
7.3
.778
10.4
.932
19.3
.963
7.4
.789
10.5
.933
19.8
.964
7.5
.799
10.7
.934
20.3
.965
7.6
.809
10.9
.935
20.9
.966
7.7
.818
11. 1
.936
21.7
.967
7.8
.826
11.3
.937
22.5
. yoo
7.9
.834
11.5
.938
23.3
.969
8.0
.842
11.7
.939
23.9
.969
1 To 4.0- inch top
in trees
6.0 inches d.b.h.
and larger.
Factor for an unlisted diameter equals factor for next smaller listed
diameter. For example, factor for 15.6 inches is .955.^
- 4 -
Table 4. — Factors for conversion of stand volumes in total cubic feet to board feet
Scribner rule per acre,1 Black Hills ponderosa pine
Average
Ratio of
Average
Ratio of
Average
Ratio of
Average
Ratio of
stand
board feet
stand
board feet
stand
board feet
stand
board feet
diameter
to
diameter
to
diameter
to
diameter
to
( Inches )
CllblC f £6t
(Inches)
cubic feet
( Inches )
cubic f S£t
( Inche s )
cubic feet
8.0
0. 78
11.9
3.
49
15. 8
4.71
19. 7
5 .
42
8.1
.85
12
0
3.56
15. 9
4.73
19.8
5 .
44
8.2
. 92
12
1
3.
61
16.0
4.76
19. 9
5.
45
8.3
.99
12
2
3.
65
16.1
4.78
20.0
5.
46
8.4
1.06
12
3
3.69
16.2
4.81
20.1
5.
47
8 . 5
1. 13
12
4
3.
73
16.3
4.83
20. 2
5.
48
8. 6
1.20
12.5
3.
77
16.4
4.86
20. 3
5.
50
8.7
1. 27
12
6
3.
80
16. 5
4.88
20.4
5.
51
8.8
1.34
12
7
3.
84
16.6
4.90
20.5
5.52
8.9
1.41
12
8
3.
88
16.7
4.92
20.6
5.
53
9 . 0
1 . 48
12
9
3.
91
16.8
4. 94
20.7
5
54
9 . 1
1 . 55
13
0
3.
95
16. 9
4.96
20.8
5 ,
56
9. 2
1 . 62
13
1
3.
98
17 . 0
4. 98
20. 9
5
57
9.3
1.68
13
2
4.
02
17.1
5.00
21.0
5.
58
9.4
1.75
13
3
4.
05
17.2
5.02
21.1
5.
59
9.5
1.82
13
4
4.
08
17.3
5.04
21.2
5.
60
9.6
1. 89
13
5
4.
11
17.4
5.06
21.3
5.
61
9. 7
1. 96
13
6
4.
14
17.5
5.08
21.4
5
62
9.8
2.03
13
7
4.
17
17.6
5.10
21.5
5
63
9.9
2.10
13
8
4
20
17.7
5.12
21.6
5
64
10. 0
2.17
13
9
4.
23
17.8
5.13
21.7
5
65
10 . 1
2 . 24
14
0
4.
25
17.9
5.15
21.8
5
66
10 . 2
2 . 31
14
1
4.
28
18.0
5.17
21.9
5
67
10.3
2.38
14
2
4
31
18.1
5.19
22.0
5
68
10.4
2.45
14
3
4.
34
18.2
5.21
22.1
5
69
10 . 5
9 ^9
Z . JZ
14
4
4.
37
18. 3
5.22
22.2
5
70
10. 6
Z . j y
14
5
4
39
18.4
5.24
22.3
5
71
10. 7
9 A1^
Z . DJ
14
6
4
42
18.5
5.26
22.4
5
72
10.8
2.72
14
7
4
44
18.6
5.27
22.5
5
73
10.9
2.79
14
8
4
47
18.7
5.29
22. 6
5
74
XX . U
Z . oo
14
9
4
49
18.8
5.30
22.7
5
75
ii l
XX . X
9
z . y j
15
0
4
52
18. 9
5.32
22.8
5
76
11.2
3.00
15
1
4
54
19.0
5.33
22.9
5
77
11.3
3.07
15
2
4
56
19 . 1
j . J j
9 ^ n
Z J . u
5
78
11.4
3.14
15
3
4
58
19.2
5.36
23.1
5
79
11.5
3.21
15
4
4
61
19.3
5.37
23.2
5
80
11.6
3.28
15
5
4
64
19.4
5.39
23.3
5
81
11.7
3.35
15
6
4
66
19.5
5.40
23.4
5
82
11.8
3.42
15
7
4
68
19.6
5.41
23.5
5
83
To 8-inch top in trees 10.0 inches d.b.h. and larger.
- 5 -
is 10 inches or larger. Other stocking levels
are named the same way. For example,
stocking level 100 means that reserve basal
area will be 100 square feet when d.b.h. after
cutting is 10 inches or larger. Basal areas
for level 100 and for diameters smaller than
10 inches are obtained by multiplying each
basal area of level 80 by the amount 100/80.
Values for other stocking levels, perhaps from
50 to 160, are computed similarly.
Periodic increases in d.b.h. due to tree
growth are estimated by the equation for DBHO
in the loop headed DO 23. The equation used
in the example is for Black Hills ponderosa
pine, and for a projection period of 10 years.
Equations for other species orprojectionperiods
may be inserted as desired. Intervals between
intermediate cuttings are one or more projection
periods long.
It is often desirable to make simulations
more realistic through introduction of variability
in values estimated by equations or contained
in tabulations. For example, repeated compu-
tations of DBHO without change in values of
the independent variables will always give the
same numerical result. In reality, actual and
estimated values frequently differ. A way of
providing variability in estimates of DBHO is
contained in the program segment between
statements 100 and 110. Similar statements
could be written for other variables.
Variability is obtained in three steps; (1)
generation of a pseudorandom number, (2) use
of this number as an independent variable to
compute the value of a residual (range: -0.3
to +0.3 inch), and (3) addition of the residual
to the computed value of DBHO. The pseudo-
random number generator, statement 100, is of
the form:
X. = AXi_1 + C (modulo M)
(Greenberger 1961). Values of all elements of
the generator are specified except for Xj_j ,
which is read in as variable GNTR. The state-
ment to compute RES is an empirical distribu-
tion function obtained by fitting a polynomial
to the normally distributed residuals of the
DBHO equation (Evans et al. 1967). An approx-
imation to the normal distribution function
may also be used (Burr 1967).
Rethinnings increase d.b.h. an average of
0.4 inch when stand densities approximate
growing stock levels most likely to be goals
of timber management.
The program can be adapted readily for
simulations with species other than the one
used for the test problem. Replacement of
tables 1 to 4, several statements of YIELD, and
species designations in two table headings are
all the substitutions necessary if the same
projection procedure is used. Statements to
be replaced are the equations for TOTO, TOTT,
DBHP, PDBHE, and DBHO. Projections based
on other procedures, such as direct estimation
of volume, can be written as a new subroutine
YIELD. It is necessary only to transfer values
of CFMO(I), BDFO(I), CFMC(I), and BDFC(I)
to the next subroutine, and to make appropriate
changes in the READ statements in INPUT1.
Subroutine ANVOL
Subroutine ANVOL is called once each test
to compute volume per acre for each year from
initial thinning to maximum stand age. Volumes
in cubic feet and board feet are computed by
linear interpolation and printed in a composite
table. Stand ages cannot exceed 179 years
unless dimensions of the 180-location arrays of
acres and annual volumes are increased.
The last few statements of the routine
expand the array of volumes cut to assign the
volume of each intermediate cut to each of the
years before the next cut is made. These
amounts are added to potential postthinning
volumes in subroutine YEARS to compute
volumes per acre of any stands not given
intermediate cutting because they are older than
minimum age for harvest.
Subroutine INPUT2
INPUT2 is called once each game to enter
numerical values of variables that may differ
for each game of a test. Descriptive data
include area of the working circle, distribution
of area by age classes, nonstocked area, and
number of acres to be planted annually.
Various costs and the rate at which they
change from year to year are entered into
computer memory. One to ten combinations
of price limit, allowable cut, and minimum
cutting age are read in for the determination
of annual cut described in the section on the
data deck.
Values read by INPUT2 that do not appear
in other tables are printed as a record of game
conditions on pages headed "alternatives for
this game."
- 6 -
Subroutine AREAS
AREAS is called once each game to compute
volumes and area distributions at the end of
the year before simulation begins. Acres in
each 1-year age class are expanded to obtain
a record of the age of each individual acre.
Total area (LAND) cannot exceed 1,000 acres
unless the dimension of ACAGE(I) is increased.
Age of the oldest acre cannot exceed 179 years
unless dimension changes listed in the descrip-
tion of AN VOL are made. Acres are totaled
by 10-year age classes. A table of initial dis-
tribution of acres by 1-year and 10-year age
classes is printed.
Growing stock volume is totaled in board
feet and in cubic feet. Volume of an acre will
be added to the total of only one of the two
volume units. The unit will be board feet if
the board-foot volume equals or exceeds the
value of variable BFMRCH read by INPUT1.
No volume will be computed for the acre if
stand age is less than the specified minimum
(AGMRCH).
Values of volume and money variables
that are not zero at start of simulation are
computed. These values are then assigned
space in one of two 2-dimensional arrays to
preserve them for printing by OUTPUT2.
Subroutine YEARS
YEARS is called each year of each game to
simulate the changes in volume and value
produced by tree growth, cutting, and other
events. Changes and their order of occurrence
are indicated by the comment statements of
the source program. Several items are described
in more detail below.
A specified area (IPLNT) is planted each
year if nonstocked acres exist. Nonstocked
acres are those deforested by fire or other
catastrophe, and do not include harvested acres
that will restock in the allotted time. Some
or all harvested acres could be added to non-
stocked area to simulate delays or failures in
natural regeneration.
Age of each acre destroyed and added to
nonstocked area is selected at random with a
pseudorandom number generator of the form:
Xi = AXi_1 + C (modulo M)
(Greenberger 1961). All values are preset except
for X j _ j which is designated as variable
ANUL. The generator has a periodicity of
128. Values of ANUL from 0 to 127 may be
selected at random to vary the pattern of loss.
Annual harvest in acres equals the constant
or variable allowable cut less any losses of
whole acres. Volume and value of shelterwood
or seed trees, if any, are credited to the year of
final cut FINL years after the regeneration
cycle starts. The volume may increase or
decrease during the regeneration period, and
may be left unharvested. Desired results are
obtained by entry of appropriate values for
GROW and SHELT.
Subroutine YEARS was written to contain
a series of dynamic events useful for many
species and forest regions. The computations
can easily be added to or modified to meet
local needs or to test special alternatives. Un-
wanted alternatives in the program need not be
removed. They can be bypassed by entry of
appropriate values for variables not needed.
Subroutine OUTPUT1
Subroutine OUTPUT1 is called after the
first year of each game and at the end of each
decade. Numbers of acres by 1-year and 10-
year stand age classes are printed. The tables
correspond to that printed by AREAS just
prior to start of a game.
Subroutine ANUAL
ANUAL is called every year of each game
to compute 40 volume, area, or money totals
and to store them for later use. Each total
is stored in one of two 2-dimensional arrays.
The first dimension identifies the variable, the
second the year of a game to which the value
applies. Numerical value of each year subscript
is year plus one so year zero of a game can
be included in the array. Array values are
used in all subsequent subroutines of MAN AGD.
Subroutine OUTPUT2
OUTPUT2 is called at the end of each game
to print the results of each year of the game.
Array values computed and stored by ANUAL
are printed in 40 numbered columns that extend
across four pages of standard Z-fold paper.
Entries under column headings are printed at
- 7 -
the rate of 40 lines, or years, per page. Column
headings on the pages produced by the test
problem (appendix 2) and the variable lists in
the source program of ANUAL (appendix 1)
identify the variables reported.
Subroutine WORTH
Subroutine WORTH is called at the end of
each game to discount all costs incurred and
all income received. Value of the growing
stock at the end of the simulation period is
discounted to beginning of the period. The
program discounts each future value at each
of 20 compound interest rates. Rates range
from 1.0 to 10.5 percent at intervals of 0.5
percent. The limits and interval can be changed
by replacement of statements for CRATE(I)
and CRATE(K) near the beginning of the
subroutine.
WORTH prints a table that gives the present
value of each of the following for each discount
rate: (1) future growing stock, (2) all incomes,
(3) sum of items 1 and 2, (4) all costs, and
(5) item 3 minus the sum of item 4 and the
value of the growing stock at beginning of the
game. Net discounted revenues (present
worths, item 5) may be plotted over discount
rates to determine the internal rate of return
applicable to the duration and conditions of the
game.
Subroutine SUMRY
Subroutine SUMRY may be called at the
end of each test to summarize results of the
games of the test. If this option is used,
SUMRY is also called at the end of each game
to store specified volume and money values in
a 3-dimensional array. Values stored correspond
to the columns of OUTPUT2 that have the
column numbers entered as KOL(I) by INPUT1.
Any of the 40 numbered columns of OUTPUT2
(appendix 2) may be reproduced. Not more
than six columns may be summarized during one
test unless the dimensions of variables KOL(I)
and SUMM(I,J,K) are increased. Results of as
many as 10 games may be summarized at one
time.
Results of the games of a test are printed
together, with a separate page for each variable
selected in advance.
Data Deck
Twenty-two types of punch cards, listed
below, are used to enter initial values of vari-
ables into computer memory. Most cards are
not optional and must be included in the data
deck so READ statements will be executed
properly. Four types are optional (6, 9, 10, 17)
and are omitted from the data deck if the
options are not to be exercised.
Data cards are read by four routines in the
order in which the types are numbered. Two
types are read once by the main program: (1)
card type 1 or BATCH name, and (2) card type
2 or the number of tests to be performed in a
batch. These identify the job and control the
number of times the rest of the main routine
is repeated.
Card types 3 to 10, inclusive, are read by
INPUT1. One card of each type except types
6, 9, and 10 must be read once each test. Card
types 6 (1 card), 9 (15 cards), and 10 (15 cards)
are omitted from the data deck if their options
are not to be used. The READ statement for
card type 6, column numbers of items to be
summarized by SUMRY, is bypassed when zero
is punched for NKOLS on card type 4. Non-
zero stumpage prices (BDPRI and/ or CFPRI)
on card type 8 cause the corresponding READ
statements for variable prices of card types
9 or 10 to be skipped.
Subroutine YIELD reads four card types
once each test. Types 11 to 14, inclusive,
contain the values of tables 1 to 4 of this
publication, or equivalent information applicable
to other species or utilization standards. A
type consists of three to eight punch cards.
Change in the lowest site index of card type
12 (table 2) or in the minimum diameter of
types 13 or 14 (tables 3 and 4) will require
changes in the statements that compute array
subscripts and probably in dimensions of the
arrays.
Card types 15 to 22, inclusive, are read by
INPUT2 once each game. Each type consists
of one card except for optional type 17, which
requires 10 punch cards. Statements that refer
to card type 17, variable area by age classes,
are bypassed when a non-zero value is punched
for KAREA on card type 16.
Card types 18, 19, and 20 contain values for
the price control procedure of Gould and O'Regan
(1965). The number of acres harvested annually
and the minimum cutting age can be made to
vary with the current stumpage price of saw
- 8 -
logs. For example, in the second game of the
test problem (appendix 2), 5 acres will be cut
if price per thousand board feet does not
exceed $12.00. Seven acres will be cut if stump-
age price is $12.01 to $15.00, and 10 acres will be
cut if price exceeds $15.00 but is less than
$99.00. Sequence of harvest is from oldest
acre to youngest, so full allowable cut will be
taken only if sufficient acres above minimum
cutting age are available. If price control is
not wanted, entries for allowable cut in columns
1 to 4 of card type 19 and for cutting age in
columns 1 to 8 of card type 20 are the desired
constant limits. A critical price greater than
the largest possible price (for example, $99.00)
is entered in columns 1 to 8 of card type 18.
Order and Contents of the Data Deck
Card Read Frequency No. of Variable
type Optional by read cards name Columns Format Description of variable
1 NO Main Batch 1 BATCH (I) 1-24 3A8 Descriptive name to identify output of
one pass through the computer.
2 NO Main Batch 1 NTSTS 1-4 14 Number of tests in the batch, each with
a yield table.
3 NO INPUTl Test 1 DESCR(I) 1-40 5A8 Phrase to describe conditions of one
test; to identify output.
4 NO INPUTl Test 1 NGAME 1-4 14 Number of trials (games) to be operated
in one test.
NKOLS 5-8 14 Number of columns of OUT PUT 2 to be
printed by SUMRY.
NOYRS 9-12 14 Number of years simulated in each game.
Can be up to 150, but will usually be
less .
INPUTl Test 1 AGEO 1-8 F8.3 Stand age at time of initial thinning.
First age given in the yield table.
SITE 9-16 F8.3 Site index. Base age and crown class
same as used to derive growth equations.
DENO 17-24 F8.3 Number of trees per acre at age AGEO.
DBHO 25-32 F8.3 Average diameter breast high of the
stand at age AGEO.
ROTA 33-40 F8.3 Maximum age in the yield table; 1 year
more than the maximum age expected dur-
ing simulations. Cannot exceed 180
years .
PRET 41-48 F8.3 Estimated percentage of the number of
trees to be retained in initial thinning
at age AGEO. Enter as a percent, e.g.,
35.0.
DLEV 49-56 F8.3 Density level for intermediate cuts
after initial thinning. Based on table
1 of this publication and procedure
described in YIELD.
CYCL 57-64 F8.3 Interval between intermediate cuts.
Equal to or a multiple of RINT.
RXNT 65-72 F8.3 Number of years for which a growth pro-
jection is made by the equation in YIELD
THIN 73-80 F8.3 Density level after initial thinning at
age AGEO. Based on table 1 and proce-
dure described in YIELD. May equal DLEV.
- 9 -
Card Read Frequency No. of Variable
type Optional by read cards name Columns Format
Description of variable
6 YES INPUT 1 Test
7 NO INPUT1 Test
NO INPUTl Test
K0L(1)
1-4
14
K0L(2)
5-8
14
K0L(3)
9-12
14
K0L(4)
13-16
14
K0L(5)
17-20
14
K0L(6)
21-24
14
AGMRCH
1-8
F8,
BFMRCH
9-16
F8,
BFSALV
17-24
F8,
COMCU
25-32
F8,
COMBF
33-40
F8.
BFPCT
41-48
F8,
CFPCT
49-56
F8.
GNTR
57-64
F8.
BDPRI
1-8
F8
CFPRI
9-16
F8
Numbers of the columns of 0UTPUT2 to be
printed by SUMRY. Maximum number of
columns is 6 and must agree with NKOLS
of card type 4. Column numbers are from
1 to 40, as given in the column headings
of the printout of 0UTPUT2 of the test
problem.
Minimum stand age for an acre to be in-
cluded in growing stock volume.
Minimum volume in M bd. ft. for an acre
to be included in board-foot growing
stock volume.
Minimum volume per acre in M bd. ft. for
commercial salvage after fire, wind, or
other loss.
Minimum cut per acre in merchantable
cubic feet (table 3) for a cut to be of
positive commercial value.
Minimum cut per acre in M bd. ft.
(table 4) for a cut to be of positive
commercial value.
Ratio, as a decimal, of board-foot
stumpage values of thinnings to board-
foot stumpage values of harvests.
Ratio, as a decimal, of cubic-foot
stumpage values of thinnings to cubic-
foot stumpage values of harvests.
Any number between 0 and 1023 used to
generate random element of the increase
from DBHT to DBHO. Enter number larger
than 1024 to bypass this step.
Stumpage price per M bd. ft. of final
harvest if price is constant for all
years of a game. Enter zero if variable
prices will be entered with card type 10.
Stumpage price per 100 cubic feet of
final harvest if price is constant for
all years of a game. Enter zero if vari-
able prices will be entered with card
type 9.
9 YES INPUTl Test
PRICF(I) 1-80 10F8.3 Stumpage price per 100 cubic feet of
harvest for each of 150 years. Used
when CFPRI equals zero.
10 YES INPUTl Test
15 PRIBD(I) 1-80 10F8.3 Stumpage price per M bd. ft. of harvest
for each of 150 years. Used when BDPRI
equals zero.
11 NO YIELD Test
TABLl(K) 1-63 21F3.1 Basal area after thinning in relation to
stand diameter. Values of table 1 copied
on punch cards. Used with DLEV and THIN.
12 NO YIELD Test
TABL2 (K,L)
1-75
25F3.1 Tree heights by age and site index class.
Values of table 2 copied on punch cards.
- 10 -
Card Read Frequency No. of Variable
type Optional by read cards name Columns Format Description of variable
13 NO YIELD Test 8 TABL3 (K) 1-72 24F3.3 Factors for conversion of total cubic
feet to merchantable cubic feet. Values
of table 3 copied on punch cards.
1A NO YIELD Test 6 TABLA(K) 1-78 26F3.2 Factors for conversion of total cubic
feet to board feet. Values of table 4
copied on punch cards.
15 NO INPUT 2 Game 1 GMNAM(I) 1-24 3A8 Descriptive name to identify each game
of a test.
16 NO INPUT 2 Game 1 LAND 1-4 14 Total acres in simulated working circle.
Maximum is 1,000 acres.
MOLD 5-8 14 Age of oldest stand in the working cir-
cle at start of a game. Maximum is 179
years .
NONSTK 9-12 I 4 Number of acres non-stocked at start of
a game. Does not include acres harvest-
ed the year before simulation begins if
regeneration will take place in the
allotted time.
KAREA 13-16 14 Number of acres in each 1-year age class
when there is equal area in each class
except for NONSTK.
IPLNT 17-20 14 Number of acres of NONSTK regenerated
annually by direct seeding or planting
at a cost of CPLT per acre.
17 YES INPUT2 Game 10 IACRE(I) 1-72 1814 Acres in each 1-year age class from 0 to
not more than 179. Use if constant area
KAREA is not wanted. Include NONSTK in
IACRE(l) as well as on card type 16.
18 NO INPUT2 Game 1 PRIDIV(I) 1-80 10F8.3 Limiting prices used to determine annual
cut in acres and minimum cutting age.
19 NO INPUT2 Game 1 MALCUT(I) 1-40 1014 Allowable annual cut in acres. May
vary with PRIDIV.
20 NO INPUT2 Game 1 FMRCHD(I) 1-80 10F8.3 Minimum cutting age. May vary with
PRIDIV.
21 NO INPUT 2 Game 1 SHELT 1-8 F8.3 Volume of shelterwood in M bd. ft. Enter
zero if shelterwood or seed trees are
not retained.
RATE 9-16 F8.3 Rate of annual increase in costs. Enter
zero if constant costs are desired.
Otherwise, enter percentage as a decimal.
CPLT 17-24 F8.3 Cost of regenerating 1 acre by seeding
or planting.
CTHN 25-32 F8.3 Cost per acre of precommercial thinning
with stand conditions as specified for
the simulation.
CLOSS 33-40 F8.3 Cost per acre of cleanup after loss due
to fire, wind, etc., when volume that
can be salvaged is less than BFSALV.
-11-
Card Read Frequency No. of Variable
type Optional by read cards name Columns Format
Description of variable
41-48 F8.3
Total per acre for 1 year of the annual
costs that can be assessed by area.
F8.3
Total of the costs that can be assessed
against each 100 cubic feet harvested.
57-64 F8.3
Total of the costs that can be assessed
against each M bd. ft. harvested.
GROW 65-72 F8.3 Percentage increase or decrease in
shelterwood volume during the regenera-
tion period. Enter as a decimal.
Enter -1.0 if shelterwood or seed trees
will not be harvested.
FINL 73-80 F8.3 Number of years between harvest cut and
removal of shelterwood or seed trees.
Enter zero if not to be removed.
22 NO INPUT2 Game 1 DEFOR 1-8 F8.5 Percentage, as a decimal, of the area of
forest lost annually to fire, wind, etc.
ANUL 9-16 F8.5 Any number between 0 and 127 used to
begin generation of pseudorandom numbers
that represent ages ot stands lost to
fire or other agency.
Description of Test Problem
The test problem that follows, (detailed in
appendix 2) demonstrates most computations
possible and the printed results obtained. It
may be used to verify accuracy of punching of
source decks and compatibility of the program
with locally available compilers. Growth projec-
tions and volume conversions are based on
relationships applicable to Black Hills ponderosa
pine (Myers 1966). Data for costs, stump age
prices, and other items are hypothetical. Results
of the simulations are therefore examples only,
and do not apply to any real forest area.
Assume an area of 915 acres of managed
stands that range from 0 (just harvested) to
129 years old. There are 7 acres in each 1-year
age class, plus 5 acres of old burn and wind-
throw that must be seeded or planted. Annual
losses to fire, wind, and other agencies average
0.04 percent of the forested area. Site index
of all acres is 60 feet.
Stands will be regenerated by two-cut
shelterwood, and will be thinned several times
during a rotation. Shelterwood volume will
average 4,000 board feet per acre, and will
increase an average of 3 percent of initial
volume each of the 10 years before the final
cut. Stands will be thinned for the first time
when they are 30 years old. At this age, there
will be 1,000 trees per acre that average 4.5
inches in diameter. Initial thinning will be to
level 120, or 120/80 times the basal areas in
table 1. Rethinnings at 20-year intervals will
be to level 100, or 100/80 times tabulated basal
areas. Stands are not expected to ever become
150 years old.
Potential prices of two products have been
estimated for each of the next 30 years. The
stumpage price of 100 cubic feet of roundwood
from mature trees or from thinnings is expected
to be $2.50 throughout the period. Price of a
thousand board feet of saw logs is expected to
vary annually, as shown in column 28 of the
printout of annual results (appendix 2). Saw
logs from thinnings will sell for 85 percent of
the price of logs from harvest cuts. A
minimum commercial cut of saw logs will be
3,000 board feet per acre, except that 1,500 board
feet may be salvaged from an acre after fire
or other catastrophe. Minimum commercial
cut of roundwood will be 400 cubic feet.
- 12 -
Current value of the growing stock will be
computed only for stands at least 40 years old.
Value will be computed for cubic volume for
acres with less than 1,500 board feet. Other-
wise, board-foot volumes will be used.
Present costs of various operations are as
follows:
Costs per acre-
Seeding— $30
Precommercial thinning— $25
Cleanup where salvage is not pos-
sible—$25
Annual costs— $0.20
Costs assessed against volume sold—
Per 100 cubic feet— $0.05
Per thousand board feet— $1.56
These costs are expected to increase at a rate
of 1 percent annually. Resources are available
to seed 1 acre each year.
Two possible means of setting the allowable
annual cut are to be tested. One alternative
is to harvest 7 acres annually, regardless of
price fluctuations. Stands less than 130 years
old will not be cut. A second possibility is to
harvest: (1) 5 acres if stumpage price per
thousand board feet is $12 or less, (2) 7 acres
if the price is $12.01 to $15, and (3) 10 acres
if price exceeds $15 per thousand. Stands less
than 130 years old will not be cut except that
the minimum age will be 120 years when
stumpage price exceeds $15. Periodic production
in board feet and net income will be compared.
Values needed to obtain present worths will
be computed.
Data cards to enter these values into compu-
ter memory must contain the alphameric
characters given in the following list. Decimal
points are shown for numbers in F-format to
indicate the way in which percentages and
money values are entered. Spaces between
numbers do not correspond to the blank columns
of the data cards. Card types 9 and 17 are not
included in the data deck because the options
that require them will be bypassed.
Test conditions and results of the simula-
tions are printed on seven types of pages
(appendix 2). The first two types, (1) a yield
table and (2) tables of volumes per acre for
each year of stand age, appear once because
one test was run. Four types of pages are
printed for each of the two games. The
seventh type of page appears once at the end
of the printout to summarize specified results
of the two games.
The two sheets of "alternatives for this
game" show the values used in the simulations,
including the different allowable cuts and cutting
ages tested.
Distributions of acres by age classes (page
type four) appear on two sets of pages, one set
for each game. Pages for year zero show 7
acres in each of 129 1-year age classes. Age
class zero has an additional 5 acres of non-
stocked area. Acreages are the same for both
games, because initial distributions were the
same. Similar pages are printed at the end of
the first year of each game and at the end of
each decade. For brevity, only the page printed
after the thirtieth year of each game is repro-
duced in appendix 2. After 30 years of simula-
tion, losses and direct seeding have modified
the pattern of 7-acre units. In addition, area
distribution of the second game has been
changed by the variable annual cuts.
The fifth type of page is a set of four pages
for each game. Values in many of the 40
numbered columns differ between games.
Board-foot volumes are unequal because of
variations in annual cuts of mature timber
during the second game. This caused money
values to differ from those reported for the
first game.
A page of discounted money values, the
sixth type of page, is printed for each game.
Rate of return was about the same for both
games. Both operations were profitable. In
addition, the forest would probably be in good
condition to produce other products, especially
recreation.
Last, specified values from each game were
printed together for convenience in interpreta-
tion of results. Total volume in board feet of
all cuts plus growing stock (column 10) was
higher after 30 years where equal areas were
cut each year. Differences were never great;
variable annual cuts produced the larger
volume after 20 years. Total net worth (column
40) was greater where annual cuts varied with
price, except for several of the earlier years.
It must be emphasized that results of these
or other simulations depend on: (1) duration of
the games, (2) values entered for the various
variables, (3) assumptions made, and (4) degree
to which the model represents reality.
The above information, additional data, and
knowledge of local conditions would help the
forest manager decide how he might best con-
duct his business. Money yields would en-
courage the manager to vary annual cuts in
- 13 -
Data Deck for Test Problem
Card
type
Contents of Cards
1
2
3
4
5
6
7
8
10
11
12
13
14
15
16
18
19
20
21
22
15
16
18
19
20
21
22
TEST PROBLEM
1
MANAGED, THINNED AGE 30
2 2 30
30. 60. 1000. 4.5
10 40
40. 1.5 1.5 400.
0.0 2.50
(Column 28 of printout by 0UTPUT2)
(Table 1 of this publication)
(Table 2 of this publication)
(Table 3 of this publication)
(Table 4 of this publication)
EQUAL AREAS CUT ANNUALLY
915
99.
7
130.
4.0
129
0
0
0
01
.0004 21.
VARY CUT WITH PRICE
915
12.
5
130.
4.0
129
15.
7
130.
.01
5
99.
10
120.
30.
7
0
0
0
25.
150.
3.0
1
0
0
0
25.
1
0
0
0
25.
40.
.85
.20
.20
100.
.05
.05
20.
1.0 2222.
0
0
0
1.56
0
0
0
1.56
10.
.30
.30
120.
0
0
0
10.0
0
0
0
10.0
0004 21.
- 14 -
response to changes in stumpage price. Highly
variable annual cuts and equally variable net
incomes would suggest that additional simula-
tions be run to test other alternatives. Cost
of computer time need not restrict the manager
in his search for information. The test problem
was compiled and run in 26 seconds of central
processor time and 9 seconds of input-output
time.
Literature Cited
Burr, Irving W.
1967. A useful approximation to the normal
distribution function, with application to
simulation. Technometrics 9: 647-651.
Chappelle, Daniel E.
1966. Economic model building and compu-
ters in forestry research. J. Forest.
64: 329-333.
Chorafas, Dimitris N.
1965. Systems and simulation. 503 pp.,
illus. New York: Academic Press.
Churchman, C. West, Ackoff, Russel L., and
Arnoff, E. Leonard.
1957. Introduction to operations research.
645 pp. New York: John Wiley and
Sons, Inc.
Clutter, Jerome L., and Bamping, James H.
1966. Computer simulation of an industrial
forestry enterprise. Soc. Amer. Forest.
Proc. 1965: 180-185.
Dane, C. W.
1966. Still more operations research. Part
III. Systems and simulation. Forest
Industries 93(11): 36-38.
Evans, George W. Ill, Wallace, Graham F., and
Sutherland, Georgia L.
1967. Simulation using digital computers.
198 pp., illus. Englewood Cliffs, New
Jersey: Prentice-Hall, Inc.
Gould, Ernest M. Jr., and O'Regan, William G.
1965. Simulation, a step toward better forest
planning. Harvard Forest Pap. 13, 86
pp., illus.
Greenberger, M.
1961. Notes on a new pseudorandom number
generator. Ass. Comput. Mach. J. 8:
163-167.
Myers, Clifford A.
1966. Yield tables for managed stands with
special reference to the Black Hills.*
U. S. Forest Serv. Res. Pap. RM-21,
20 pp., illus. Rocky Mountain Forest and
Range Exp. Sta., Ft. Collins, Colo.
and Godsey, Gary L.
1968. Rapid computation of yield tables for
managed, even-aged timber stands.*
U. S. Forest Serv. Res. Pap. RM-43,
16 pp., illus. Rocky Mountain Forest
and Range Exp. Sta., Ft. Collins Colo.
O'Regan, William G., Arvanitis, Lucas, and
Gould, Ernest M. Jr.
1966. Systems simulation, and forest man-
agement. Soc. Amer. Forest. Proc.
1965: 194-198.
Thompson, Emmett F.
1966. Traditional forest regulation model:
An economic critique. J. Forest. 64:
750-752.
♦Address requests for copies to the originating
office.
- 15
APPENDIX I
Program MANAGD
DEFINITIONS OF VARIABLES
ACAGEII) = ONE ACRE, WITH AGE OF NUMERICAL VALUE AND IDENTIFIED BY
SUBSCRIPT
ACCST = ANNUAL COST PER ACRE
AGEO = INITIAL AGE IN YIELD TABLE
AGMRCH = MINIMUM AGE FOR STAND TO BE INCLUDED IN GROW I NG STOCK
ANBDF(I) = M BD. FT. PER ACRE AT END OF EACH YEAR
ANCUV(I) = CU.FT. STANDING PER ACRE AT END OF EACH YEAR
ANNE T = ANNUAL NET INCOME
ANUL = NUMBER BETWEEN 0 AND 127 USED TO START- GENERATION OF
PSEUDORANDOM NUMBERS
BASC = BASAL AREA REMOVED PER ACRE
BASO = BASAL AREA PER ACRE BEFORE THINNING
BAST = BASAL AREA PER ACRE AFTER THINNING
BATCH ( I ) = JOB NAME
BDFC(I) = M BD. FT. REMOVED PER ACRE
BDFO(I) = M BD. FT. PER ACRE BEFORE THINNING
BDFT(I) = M BD. FT. PER ACRE AFTER THINNING
BOPRI = CONSTANT STUMPAGE PRICE PER M BD. FT.
BFCST = COSTS PER M BD. FT. HARVESTED
BFMRCH = MINIMUM VOLUME TO BE INCLUDED IN BD. FT . GROWING STOCK
BFPCT = PCT. TO CONVERT BD.FT. PRICE FOR THINNINGS
BFSALV = MINIMUM BD.FT. FOR COMMERCIAL SALVAGE
CFMC(I) = MERCHANTABLE CU.FT. REMOVED PER ACRE
= MERCHANTABLE CU.FT. PER ACRE BEFORE THINNING
= MERCHANTABLE CU.FT. PER ACRE AFTER THINNING
PCT. TO CONVERT CU.FT. PRICE FOR THINNINGS
CONSTANT STUMPAGE PRICE PER 100 CU.FT.
COST OF CLEANUP OF VOLUME NOT SALVAGED
MINIMUM COMMERCIAL CUT IN BD. FT.
MINIMUM COMMERCIAL CUT IN CU. FT.
PLANTING COST PER ACRE
CFMO ( I )
CFMT ( I )
CFPCT =
CFPRI =
CLOSS =
COMBF =
COMCU =
CPLT =
CRATE! I) = INTEREST RATES FOR DISCOUNTING
CSTAC = ANNUAL COSTS BASED ON AREA
CSTVL = ANNUAL COSTS FOR VOLUME HARVESTED
CTHN = COST PER ACRE OF PR ECOMME RC I AL THINNING
CUCST = COSTS PER 100 CUBIC FEET HARVESTED
CUTAGE = MINIMUM CUTTING AGE
CYCL = INTERVAL BETWEEN INTERMEDIATE CUTS
DBHO = AVERAGE STAND D.B.H. BEFORE THINNING
DBHT = AVERAGE STAND D.B.H. AFTER THINNING
DEFOR = PERCENTAGE! AS DECIMAL) OF NUMBER OF ACRES LOST ANNUALLY
DENC = TREES REMOVED PER ACRE
DENO = TREES PER ACRE BEFORE THINNING
DENT = TREES PER ACRE AFTER THINNING
OESCR(I) = DESCRIPTION OF TEST CONDITIONS
DISCI I) = DISCOUNTED VALUE OF FUTURE COSTS
DISG(I) = DISCOUNTED VALUE OF GROWING STOCK
DISI(I) = DISCOUNTED VALUE OF FUTURE INCOMES
DLEV = GROWING STOCK LEVEL FOR SECOND AND SUBSEQUENT THINNINGS
FINL = YEARS BETWEEN HARVEST AND REMOVAL OF SHELTERWOOD
FMRCHD ( I ) = MINIMUM CUTTING AGE BASED ON PRICE
GMNAM(I) = NAME OF THE GAME
GNTR = PSEUOORANOOM NUMBER GENERATOR. VALUE 0 TO 1023.
GROW = GROWTH RATE OF SHELTERWOOD
GSVALB = DOLLAR VALUE OF BD.FT. GROWING STOCK
GSVALC = DOLLAR VALUE OF CU.FT. GROWING STOCK
GVLBF = GROWING STOCK VOLUME, M BD. FT.
GVLCU = GROWING STOCK VOLUME, CU.FT.
HTSO = TREE HEIGHT BEFORE THINNING
HTST = TREE HEIGHT AFTER THINNING
IACKE(I) = ACRES OF WORKING CIRCLE IN EACH 1-YEAR AGE CLASS, AT
START OF GAME
IALCUT = NUMBER OF ACRES ALLOWABLE ANNUAL CUT
IGAME = NUMBER OF GAME
IPLNT = NUMBER OF NON-STOCKED ACRES REGENERATED ANNUALLY
ISUM(I) = TOTAL ACRES EACH 10-YR AGE CLASS
ITEST = NUMBER OF TEST
IVAR(I.J) = VARIABLES TO BE PRINTED BY 0UTPUT2
IYEAR = YEAR WITHIN RUN OF A GAME
KAREA = EQUAL AREA OF EACH 1-YEAR AGE CLASS
KOL(I) = COLUMN NUMBER ( FROM 0UTPUT2 ) PRINTED BY SUMRY
KOUNT = COUNT OF ACRES HARVESTED, PLUS ONE
LAND = TOTAL ACRES IN SIMULATED WORKING CIRCLE
LAST = NUMBER OF LAST ACRE HARVESTED
MALCUT(I) = ANNUAL ALLOWABLE CUT BASED ON PRICE
MOLD = AGE OF OLDEST ACRE IN WORKING CIRCLE AT START OF A GAME
IJACREII) = ACRES OF WORKING CIRCLE IN EACH 1-YEAR AGE CLASS,
DURING A GAME
NGAMF = NUMBER OF GAMES PER TEST
NKOLS = NUMBER OF COLUMNS OF 0UTPUT2 TO BE PRINTED BY SUMRY
NONSTK = NONSTOCKED AREA
NOYRS = NUMBER OF YEARS IN A GAME
NT ST S = NUMBER OF TESTS IN BATCH
PRFT = PERCENTAGE OF TREES RETAINED AFTER INITIAL THINNING
PREV(I) = PRESENT VALUE OF GROWING STOCK AND INCOMES
PRIBD(I) = VARIABLE STUMPAGE PRICE PER M BD. FT.
PRICF(I) = VARIABLE STUMPAGE PRICE PER 100 CU.FT.
PRIDIV(I) = PRICES USED TO SET POLICY
PWTH(I) = PRESENT WORTH
RATE = RATE OF ANNUAL INCREASE IN COSTS
RETHV = ANNUAL RETURN FROM FINAL HARVEST
RETRN = ANNUAL INCOME FROM STUMPAGE
RET TH = ANNUAL RETURN FROM THINNINGS
R I NT = NUMBER OF YEARS FOR WHICH GROWTH PROJECTION IS MADE
ROTA = LONGEST POSSIBLE ROTATION IN YIELD TABLE
SCPLT = TOTAL ANNUAL PLANTING COST
SCTHN = SUM OF PRECOMMERC I AL THINNING COSTS
SHELT = M BD. FT. RETAINED AS SEED TREES OR SHFLTERWOOD
SHUD = VOLUME OF SHELTERWOOD AT FINAL CUT
SITE = SITE INDEX
SUMM(I,J,K) = ARRAY FOR PRINTING BY SUMRY
TCOST = TOTAL ANNUAL COSTS
THIN = GROWING STOCK LEVEL FOR INITIAL THINNING
TOTC = TOTAL CUBIC FEET REMOVED PER ACRE
rOTO = TOTAL CUBIC FEET PER ACRE BEFORE THINNING
TOTT = TOTAL CUBIC FEET PER ACRE AFTER THINNING
VAR(I.J) = VARIABLES TO BE PRINTED BY 0UTPUT2
VBHV = BOARD-FOOT VOLUME FROM HARVESTS
VB TH = BOARD-FOOT VOLUME FROM THINNING
VCHV = CUBIC-FOOT VOLUME FROM HARVESTS
VCTH = CUBIC-FOOT VOLUME FROM THINNING
VLBF = VOLUME HARVESTED, M BD. FT.
VLCU = VOLUME HARVESTED, CU.FT.
YRLUS = NUMBER OF ACRES LOST ANNUALLY
COMMON ACCST, AGEO, AGMRCH , ANBDF ( 181) , ANCUVI 181) , ANNE T , BATCH ( 3 ), BDFC
1 ( 180) ,BDFO( 180) , BFCST, CFMCI 180) ,CFMO( 180), CSTAC,CSTVL, CUCST, CUTAGE
2,DBH0,DEN0,DESCR( 5) , DL EV , FMRCHD ( 10) , GMNAM I 3 ) , GSVALB , GSVALC , GVLBF ,
3GVLCU, I ACRE (180) , I ALCUT , I G AME , I SUM ( 18 ), ITEST, I VAR ( 26 , 1 50 ) , I YEAR ,
4K0L (6),LAND,LAST,MALCUT(10),M0LD,NACRE(180),NGAME,NKOLS,Nl , NONSTK,
5N(JYRS,PRET, PRIBDI 150) ,PRICF(150),PRIDIV( 10) , RETRN, ROTA, RATE, GROW,
6SITE,SUMM(6,25, 10), TCOST, SHE LT.YR LOS, FINL, CLOSS, CTHN, CPLT, CFPCT,
7VAR( 14,150) , VLBF, VLCU, ACAGE ( 1000) , CYCL, KOUNT, DEFOR, ANUL, MIX, IPLNT,
8RINT, THIN, BFMRCH, BFSALV, COMCU, COMBF, BFPCT, GNTR
READ (5,1) (BATCH(I),I=1,3)
1 FORMAT (3A8)
READ (5,2) NTSTS
2 FORMAT (14)
DO 8 ITEST=1, NTSTS
CALL INPUT1
PRINT YIELD TABLE AND COMPUTE VOLUME FOR EACH YEAR OF STAND AGE
OPERATE SYSTEM FOR DESIRED NUMBER OF GAMES
DO 8 I GAME= 1 » NG AME
CALL I NPUT2
CREATE ACRES IN EACH AGE CLASS
OPERATE SYSTEM FOR DESIRED NUMBER OF YEARS
DO 7 I YE AR= 1 , NOYRS
CALL YEARS
PRINT ACRES IN EACH AGE CLASS FOR FIRST YEAR AND END OF
EACH DECADE
IF ( I YEAR .LE. 1 ) GO TO 3
IF (UK .EQ. 10) GO TO 4
GO TO 6
3 UK = 1
GO TO 5
4 UK = 0
5 CALL 0UTPUT1
6 UK = UK + 1
CALL ANUAL
7 CONTINUE
PRINT VOLUMES AND VALUES FOR EACH YEAR
CALL 0UTPUT2
CALL WORTH
SUMMARIZE DESIRED NUMBER OF COLUMNS OF 0UTPUT2
IF (NKOLS .LE. 0) GO TO 8
CALL SUMRY
8 CONTINUE
CALL EXIT
END
SUBROUTINE INPUT1
COMMON ACCST, AGEO, AGMRCH, ANBDF ( 181 ) , ANCUVI 181 ) , ANNE T , BATCH ( 3 ) , BDFC
1 ( 180) ,BDFO( 180) , BFCST, CFMCI 180) ,CFMO( 180 ), CSTAC , CSTVL , CUCST, CUTAGE
2, DBHO, DENO, DESCR(5),DLEV,F MR CHD(10),GMNAM(3), GSVALB, GSVALC, GVLBF,
3GVLCU, I ACRE ( 180) , IALCUT, IGAME, TSUMI 18) , ITEST, IV AR (26, 150), IYEAR,
4K0L ( 6 ), LAND, LAS T, MALCUT ( 10 ), MOLD, N ACRE! 180 ),NGA ME, NKOLS.Nl, NONSTK,
5N0YR S, PRE T, PRIBDI 150) ,PRICF(150),PRIDIV(10) , RE T RN , ROT A , RATE , GROW ,
6SITE,SUMM(6,25, 10) , TCOST , SHFLT , YRLOS , F I NL , CLOS S , C THN , CPL T , CF PC T ,
7VARI 14, 150) , VLBF, VLCU, AC AGE ( 1000) ,CYCL, KOUNT, DEFOR, ANUL, MIX, IPLNT,
8R I NT, THIN, BFMRCH, BFSALV, COMCU, COMBF, 6FPCT, GNTR
SET INITIAL VALUES OF ZERO
DO 1 1=1,6
1 KOLI I ) = 0
DO 2 1=1,150
PRIBDI I I = 0.0
2 PRICF I I ) = 0.0
DO 3 1=1,181
ANBDF I I ) = 0.0
3 ANCUVI I ) = 0.0
DO 4 1=1,6
DO 4 J=l,25
DO 4 K=l, 10
4 SUMM( I , J,K) = 0.0
(EAD VALUES THAT DO NOT CHANGE WITHIN A TEST
- 16 -
= 1,5)
READ (5,5) (DESCRI I),
5 FORMAT (5A8)
READ (5,6) NGAME,NKOLS»NOYRS
6 FORMAT (2014)
READ (5,7) AGEO, SITE, DENO.DBHO, ROTA, PRE T , DL E V, C YCL , R I NT , THIN
7 FORMAT (10F8.3)
IF (NKOLS .LE. 0) GO TO 8
READ (5,6) ( KOL ( I ) , 1=1, NKOLS)
8 READ (5,7) AGMRCH , BFMRCH , B FS AL V , C OMCU , COMBF , BF PCT , C FPC T , GNTR
READ (5,7) BDPRI.CFPRI
CREATE A SERIES OF CONSTANT OR VARIABLE PRICES
IF (CFPRI .NE. 0.0) GO TO 9
READ (5,7) ( PRICFI I ) , 1=1, 150)
GO TO 11
9 DO 10 1=1,150
10 PRICFI I ) = CFPRI
11 IF (BDPRI -NE. 0.0) GO TO 12
READ (5,7) ( PRIBDI I ) , 1=1, 150)
GO TO 14
12 DO 13 1=1,150
13 PRIBDI I ) = BDPRI
14 RETURN
END
SUBROUTINE YIELD
COMMON ACCST, AGEO, AGMRCH, ANBOF I 181 ) ,ANCUV( 181) ,ANNET, BATCH! 3) ,BDFC
1(180),BDFO(180),BFCST,CFMC(180) ,CFMO( 1 80 ) , C STAC , C ST VL , C UCST , CUT AGE
2,DBH0,DEN0,0ESCR(5 ) , DLE V , FMRCHD I 1 0 ) , GMN AM I 3 ) , GS V ALB , GS V ALC , GVL BF ,
3GVLCU, I ACRE! 180) , I ALCUT , I GAME , I SUM I 18) , ITEST, IVARI26, 150) , I YEAR,
4K0L ( 6 ), LAND, LAST, MALCUTI 10), MO LD,N ACRE! 180 ),NGAME, NKOLS, NL.NONSTK,
5N0YRS, PRE T, PRIBDI 150) , PR I C F ( 1 50 ) , PR I D I V I 10), RETRN, ROTA, RATE, GROW,
6SITE,SUMM(6,25, 10 ) , TCOST , SHELT , YRLOS , F I NL , C LOS S , CTHN , C PLT , C FPC T ,
7VARI 14, 150) ,VLBF,VLCU,ACAGE( 1000) , C YCL , KOUNT , DE FOR , ANUL , M I X , I PL NT ,
8RINT, THIN, BFMRCH, BFSALV,COMCUiCOMBF,BFPCT, GNTR
DIMENSION TABL1 (81 ) ,TABL2(4,18) , T ABL3 I 192 ) , T ABL4 I 1 56 )
BDFT = 0.0
CFMT = 0.0
JBDFC = 0
J8DF0 = 0
JBDFT = 0
JCFMC = 0
JCFMO = 0
JCFMT = 0
DO 1 1=1,180
BDFCI I ) =0.0
BDFOI I ) = 0.0
CFMCI I ) = 0.0
1 CFMOI I ) = 0.0
READ HEIGHTS, DENSITIES, VOLUME CONVERSION FACTORS
READ (5,2) (TABL1IK),K=1,81)
2 FORMAT (21F3.1)
READ (5,3) ( (TABL2(K,L),K=1,4),L=1,18)
3 FORMAT (25F3.1)
READ (5,4) (TABL3(K),K=1,192)
4 FORMAT (24F3.3)
READ (5,5) (TABL4IK) ,K=1, 156)
5 FORMAT (26F3.2)
Nl = AGEO
N = AGEO
BASO = DENO » 0.0054542 » DBHO * DBHO
OBTAIN HTSO
TEN = 10.0
ISITE = (SITE/TEN - 3.0) + 0.01
I AGEO = (AGEO/TEN) + 0.01
HTSO = TABL2I ISITE, I AGEO)
TOTO = (0.4047 » BASO • HTSO) + (25.5970 » DBHO) - 191.6433
COMPUTE MERCH. CU.FT. IF DBH IS AT LEAST 5.0 INCHES
IF (DBHO .LT. 5.0) GO TO 6
I DBHO = ((DBHO - 5.0) » 10.0) + 1.01
XDBHO = TABL3I I DBHO )
CFMO(N) = TOTO • XDBHO
COMPUTE BD.FT. IF DBH IS AT LEAST 8.0 INCHES
IF (DBHO .LT. 8.0) GO TO 6
JDBHO = ((DBHO - 8.0) * 10.0) + 1.01
YDBHO = T ABL4 ( JDBHO )
BDFO(N) = TOTO * YDBHO
COMPUTE DBH AFTER INITIAL THINNING
DO 11 J=l,100
PDBHE = 0.95462»ALOG101DBHO) -
DBHE = 10.0 •* PDBHE
IDBHE = DBHE • 10.0 + 0.5
DBHE = IDBHE
DBHE = DBHE/TEN
DENE = DENO » (PRET/100.0)
BASE = (0.0054542 * DBHE * DBH
BREAK = 49.9 * THIN/80.0
IF (BASE .GT. BREAK) GO TO 7
DI3HP = (80.0/THIN) * (0.08733
GO TO 8
DBHP = (80.0/THIN) * (0.10938
IUBHP = DBHP • 10.0 + 0.5
DBHP = I DBHP
DBHP = D8HP/TEN
IF I DBHP - DBHE) 9,12,10
PRET = PRET + PRET • 0.02
GO TO 11
0.10640*ALOG10IPRET)
BASE) - 0.17858
10 PRET = PRET - PRET • 0.02
11 CONTINUE
12 DBHT = DBHE
COMPUTE VALUES AFTER INITIAL THINNING
JDBHT = ( I DBHT - 2.0) » 10.0) + 1.01
SOFT = TABL1 I JDBHT )
BAST = (THIN/80.0) « SOFT
ENTER LOOP FOR ALL REMAINING COMPUTATIONS AND PRINTOUT
DO 26 1=1,100
DENT = BAST/IO. 0054542
HTST = HTSO
TOTT = (0.4047 » BAST
COMPUTE MERCH. CU.FT.
♦ DBHT * DBHT)
* HTST) + (25.5970 * DBHT) - 191.6433
IF DBH IS AT LEAST 5.0 INCHES
IF (DBHT .LT. 5.0) GO TO 13
I DBHT = I (DBHT - 5.0) • 10.0)
XDBHT = TABL3 ( I DBHT )
CFMT = TOTT * XDBHT
1.01
COMPUTE BD.FT. IF DBH IS AT LEAST 8.0 INCHES
IF (DBHT .LT. 8.0) GO TO 13
JDBHT = ((DBHT - 8.0) * 10.0) + 1.01
YOBHT = T ABL4 ( J DBHT )
BDFT = TOTT * YDBHT
CHANGE MODE AND ROUND OFF FOR PRINTING
13 JDENO
JHTSO
JTOTO
JBASO
JCFMO
JBDFO
JBDFO
JDENT
JHTST
J TOTT
JDENC
JCFMT
CFMT
JB AST
JBASC
JTOTC
JCFMC
DENO
HTSO
TOTO + 0.5
BASO + 0.5
CFMOIN) + 0.5
( BDFO ( N ) /TEN ) + 0.5
JBDFO » 10
DENT
HTST
TOTT + 0.5
JDENO - JDENT
CFMT + 0.5
JCFMT
IF (JCFMT .GT. JCFMO) JCFMO = JCFMT
CFMOIN) = JCFMO
JBDFT = I BDFT/TEN ) + 0.5
JBDFT = JBDFT * 10
BDFT = JBDFT
BDFT = BDFT » .001
IF (JBDFT .GT. JBDFO) JBDFO = JBDFT
BDFO(N) = JBDFO
BDFO(N) = BDFO(N) • .001
BAST + 0.5
JBASO - JBAST
JTOTO - J TOTT
JCFMO - JCFMT
IF (JCFMC .LE. 0) JCFMC = 0
CFMC(N) = JCFMC
JBDFC = JBDFO - JBDFT
IF (JBDFC .LE. 0) JBDFC = 0
BDFC(N) = JBDFC
BDFC(N) = BDFC(N) * .001
IF ( I . GE . 2) GO TO 19
WRITE HEADINGS FOR YIELD TABLE
WRITE (6,14) SITE.CYCL.DLEV
14 FORMAT ( 1H1,////,28X,81HYIELDS PER ACRE OF MANAGED, EVEN-AGED STAN
IDS OF PONDEROS A PINE IN THE BLACK HILLS/1H ,35X,10HSITE INDEX, F3.0
2,1H, ,F4.0,19H-YEAR CUTTING C YC LE , 1H , , 1 5H DENSITY LEVEL ,F4.0,///)
WRITE (6,15)
15 FORMAT ( 1H0,25X,38HENTIRE STAND BEFORE AND AFTER THINNING, 28X,26HP
1ERI0DIC CUT AND MORTALITY)
WRITE (6,16)
16 FORMAT ( 1H0,9X,5HSTAND, 1 OX , 5HB AS AL , 3 X , 7H AVER AGE , 2 X , 7H A VER AGE , 3 X , 5H
1T0TAL,3X,9HMERCHANT-,3X,9HSAWTIMBER,9X,5HBASAL,4X,5HT0TAL,3X,9HMER
2CHANT-,3X,9HSAWTIMBER)
WRITE (6,17)
17 FORMAT (1H , 1 OX , 3H AGE , 4X , 5HTRE E S , 3X , 4H AR E A , 4X , 6HD . B . H . , 3X , 6HHE I GHT
1,2X,6HV0LUME,2X,11HABLE VOLUME , 4X , 6H VOLUME , 3X , 5HTRE E S , 3 X , 4H ARE A , 3 X
2, 6H VOLUME, 2X.11HABLE VOLUME , 4X , 6HV0LUME )
WRITE (6,18)
18 FORMAT ( 1H , 8X , 7H ( YE ARS ) , 3X , 3HN0 . , 3 X , 6HS0. F T . , 4X , 3H I N. , 6X , 3HFT . , 4X
1,6HCU.FT.,5X,6HCU.FT.,8X,3HMBF,5X,3HN0.,3X,6HSQ.FT.,2X,6HCU.FT.,5X
2,6HCU.FT.,8X,3HMBF)
19 WRITE (6,20) AGEO , JDENO , JBASO , DBHO , JHTSO , JTOTO , C FMO ( N ), BDFO ( N )
20 FORMAT ( 1H0 , 9X , F4. 0 , 4X , I 5 , 2X , I 4 , 5X , F 5 . 1 , 5X , I 3 , 4X , I 5 , 6X , F 5 . 0 , 6 X , F 6.
13)
IF (AGEO. GE. ROTA) GO TO 27
WRITE (6,21 ) AGEO, JDENT, JBAST , DBHT, JHTST, J TOTT, CFMT, BDFT, JDENC, J 8 A
ISC, JTOTC, CFMC IN) , BDFC I N )
21 FORMAT I 1H , 9X , F4 . 0 , 4X , I 5 , 2X , I 4, 5X , F5 . 1 , 5X , I 3 , 4X , I 5 , 6X , F5 . 0 , 6X , F6.
13.4X, I5..3X, I3.5X, I4,6X,F4.0,8X,F5.3)
COMPUTE VALUES FOR EACH PERIOD. THIN AS SPECIFIED
IK = CYCL/RINT
DO 2 3 L-ltlK
AGEO = AGEO + R I NT
N = AGEO
IF (AGEO . GT . ROTA) GO TO 27
DBHO = 1.0097»DBHT + 0.0096*SITE
I DBHO = DBHO » 10.0 + 0.5
D8H0 = I DBHO
DBHO = DBHO/TEN
1.5766*ALOG10(BAST)
- 17 -
ADD RANDOM ELEMENT TO PREDICTED DBHO.IF DESIRED
IF (GNTR .GT. 1024.0) GO TO 110
IDIV = (17.0 » GNTR + 3.01/1024.0
NGNTR = GNTR
GNTR = (17 * NGNTR + 3) - 1024 » IDIV
IF (GNTR .GT. 1000.0) GO TO 100
IF (GNTR .LT. 0.0) GO TO 100
Al = GNTR/100.0
A2
Al
0.0523
A2
IRES
IRES
0.5
0.5
ADJ/TEN
RFS = 0.9565
1*A2 - 3.3009
IRES = RES
IF (RES .LT.
IF (RES . GT.
ADJ = IRES
DBHO = DBHO
110 DENO = DENT
8AS0 = DENO * (0.0054542 »
ISITE = (SITE/TEN - 3.0) +
IAGEO = (AGEO/TEN) + 0.01
HTSO = TABL2I ISITE, IAGEO)
TOTO = 0.4047 » BASO » HTSO
COMPUTE MERCH. CU.FT. IF DBH IS AT LEAST 5.0 INCHES
IF (DBHO .LT. 5.0) GO TO 22
I DBHO = ((DBHO - 5.0) * 10.0) + 1.01
XDBHO = TABL3 ( IDRHO)
CFMO(N) = TOTO • XDBHO
COMPUTE BD.FT. IF DBH IS AT LEAST 8.0 INCHES
IF (DBHO .LT. 8.0) GO TO 22
JDBHO = ((DBHO - 8.0) * 10.0) + 1.01
YDBHO = TABL4 ( JDBHO )
BDFO(N) = TOTO » YOBHO
22 IF (L .EC. IK) GO TO 24
WRITE VALUES FOR END OF PERIOD IF THINNING NOT DUE
A2 + 0.00084
DENO
HTSO
BASO + 0.5
TOTO + 0.5
CFMO(N) + 0.?
= JCFMO
( BDFO ( N ) /TEN I
JBDFO » 10
= JBDFO
= BDFO(N) * .001
AGEO,KDENO,KBASO, DBHO, KHTSO.KTOTO, CFMO ( N ) , BDFO ( M)
0.5
KQENO =
KHTSO =
KBASO =
KTOTO =
JCFMO =
CFKOIN)
JBDFO =
JBDFO =
BDFO ( N )
BDFO(N)
WRITE (6,20)
DBHT = DBHO
BAST = BASO
23 CONTINUE
INCREASE DBH AS RESULT OF THINNING AND COMPUTE POST-THINNING
VALUES
24 DBHT = DBHO + 0.4
IF (DBHT . GE . 10.0) GO TO 25
JOBHT = ((DBHT - 2.0) * 10.0) + 1.01
SQFT = TABL 1 ( JDBHT )
BAST = (DLEV/80.0) * SQFT
GO TO 26
25 BAST = DLEV
26 CONTINUE
27 RETURN
END
SUBROUTINE ANVOL
COMMON ACCST, AGEO, AGMRCH, ANBDFI 1 8 1 1 , ANCUV ( 1 8 1 1 , ANNET , B ATCH ( 3 1 , BDFC
1 ( ieO) ,BDFO( 180) ,BFCST,CFMC( 180) ,CFMO( 1 8 0 ) , C STAC , C STVL , CUC ST , CUT AGE
2,DBHO,OENO,DESCR(5) , OLE V , FMRCHD ( 10) ,GMNAM(3) , GS VALB , GSVALC , GVLBF ,
3GVLCU, I ACRE (180) , IALCUT, I GAME, ISUMf 18),ITEST,IVAR(26,150) , I YEAR,
4K0LI6) , LAND, LAST, MALCUTI 10) .MOLD, NACRE (180) , NGAME , NKOLS , Nl , NONSTK ,
5N0YRS,PRET,PRIBD( 150),PRICF(150),PRIDIV( 10) ,RETRN, ROTA, RATE, GROW,
6SITE,SUMM(6,25, 10 ) , TCOST , SHELT , YRLOS , F I NL , C LOSS , C THN , C P LT , CF PC T ,
7VAR(14,150),VLBF,VLCU,ACAGE(1000),CYCL,K0UNT,DEF0R,ANUL,MIX,IPLNT,
8RINT,THIN,BFMRCH,BFSALV,C0MCU,C0MBF,BFPCT, GNTR
IROT = ROTA
INT = R I NT
.NVUL = ( ( IROT - Nl ) /INT) + 1
K = NVOL - 1
INTERPOLATE BETWEEN VALUES FROM YIELD TABLE
'DO 1 L=1,K
DO 1 J=1,INT
NN = J + Nl + ( L - 1 ) « INT
RJ = J - 1
N = Nl + (L - 1) » INT
ANCUV (NN) =CFMO(N)-CFMC(N)+(RJ/RINT)(
ANBDFI NN) = BDFO ( N 1 -BDFC ( N 1 + ( R J / R I NT 1 <
1 CONTINUE
(CFMO(N+INT)-CFMO(N)+CFMC(N) )
(BDFO(N+INT)-BDFO(N)+BDFC(N) )
WRITE TABLE HEADINGS
WRITE (6,2) SITE, CYCL, THIN, DLEV
2 FORMAT ( 1H1,////,41X,51HGR0WING STOCK OF MANAGED BLACK HILLS PONDE
1R0SA PINE/1H ,47X,10HSITE INDEX, F3.0,1H,,F4.0,19H-YEAR CUTTING CYC
2LE/1H ,53X, 14HDENSITY L E VEL- , F4 . 0 , IX , 3H AND , F4 . 0 )
WRITE (6,3)
3 FORMAT ( 1H0,43X,44HV0LUMES PRESENT PER ACRE AT END OF EACH YEAR)
WRITE (6,4)
4 FORMAT ( 1H ,54X,23HMERCHANTABLE CUBIC FEET/1H0,64X,4HYEAR/1H ,14X,
16HDECADE,9X,1H0,9X, 1H1.9X, 1H2 , 9X , 1H3 , 9X , 1H4 , 9X , 1H 5 , 9X , 1H6 , 9X , 1H7 , 9
2X, 1H8.9X, 1H9,// )
WRITE CUBIC FEET PER ACRE FOR EACH YEAR
WRITE (6,5) K, ( ANCUVI NN 1 ,NN=1 , 10)
i FORMAT (1H , 120, F13. 1 ,9F10. 1 )
IJ = IROT/10 - 1
DO 6 J = l , I J
NN=10*J+1
WRITE (6,5) J, ANCUV (NN) , ANCUV (NN+1) , ANCUV (NN + 2 ) , ANCU V ( NN+3 ), ANCUV (
1NN + 4) , ANCUV (NN + 5) , ANCUV ( NN + 6 ) , ANCUV ( NN + 7 ) , ANCUV (NN + 8) , ANCUV ( NN + 9 )
, CONTINUE
J = IROT/10
ANCUVI IROT+l ) = CFMO(IROT)
WRITE (6,5) J , ANCUVI IROT+l )
WRITE BD.FT. PER ACRE FOR EACH YEAR
WRITE (6,7)
' FORMAT I 1H0,///,55X,23HTH0USANDS OF BOARD FEET,//)
WRITE (6,8) K, ( ANBDF(NN) ,NN=1, 10)
I FORMAT ( 1H ,I20,F13.3,9F10.3)
DO 9 J = 1,1 J
NN = 10 » J + 1
WRITE (6,8) J,ANBDF(NN) ,ANBDF(NN + 1) ,ANB0F(NN + 2) ,ANBDF(NN + 3), ANBDFI
l.NN + 4) , ANBDF ( NN + 5 ) ,ANBDF(NN+6) ,ANBDF(NN + 7) ,ANBDF(NN + 8) , ANBDF ( NN + 9 )
I CONTINUE
J = IROT/10
ANBDFI IROT+l) = BDFOI IROT)
WRITE (6,8) J, ANBDFI IROT+l )
PROVIDE FOR ANY ACRES BEYOND CUTAGE LEFT UNTHI NNED FOR A FEW YEARS
JCYCL = CYCL
IMI = JCYCL - 1
DO 10 I=N1, IROT, JCYCL
DO 10 J=1,IMI
NX = I + J
180) GO TO 11
BDFC I I )
CFMCI I 1
IF (NX .GE.
BDFC ( NX )
10 CFMC(NX)
11 RETURN
END
SUBROUTINE INPUT2
COMMON ACCST, AGEO, AGMRCH, ANBDFI 181) , ANCUV I 181 ) , ANNET , BATCH I 3 ), BDFC
1(180), BDFOI180) ,BFCST, CFMCI 180) ,CFMO( 180 ), CST AC , C STVL , CUC ST , CUTAGE
2,DBH0,DEN0,DESCR( 5) , DLE V , F MRCHD ( 1 0 1 , GMN AM ( 3 1 , G S V A LB , GS V ALC , GVLBF ,
3GVLCU, I ACRE! 180) , IALCUT, I G AME , I SUM ( 1 8 ) , I T EST , I VAR ( 26 , 150) , I YEAR,
4K0LI6) , LAND, LAST, MALCUTI 10), MOLD, NACRE I 180) , NG AME , NKOLS , N 1 , NONSTK ,
5N0YRS, PRET, PRIBDI 150),PRICF(150),PRIDIV( 10 1 , RE TRN , ROT A , RATE , GROW ,
6SITE,SUMM(6, 25, 10), TCOST, SHELT, YRLOS, FINL,CLOSS,CTHN,CPLT,CFPCT,
7VAR(14,150) ,VLBF, VLCU, AC AGE ( 1000) , C YCL , KOUNT , DE FOR , ANUL , M I X , I PL NT ,
8RINT,THIN,BFMRCH,BFSALV,C0MCU,C0MBF, BFPCT.GNTR
SET INITIAL VALUES OF ZERO
CSTAC = 0.0
CSTVL = 0.0
IALCUT = 0
KOUNT = 1
LAST = 0
RETRN = 0.0
YRLOS = 0.0
DO 1 1=1,180
I ACRE ( I ) = 0
1 NACRE ( I ) = 0
DO 2 1=1,10
FMRCHD ( I ) = 0.0
MALCUT(I) = 0
2 PRIDIVI I ) = 0.0
DO 3 1=1,26
DO 3 J=l,150
3 IVARI I , J) = 0
DO 4 1=1,14
DO 4 J=l,150
4 VAR( I , J ) = 0.0
READ VALUES THAT DO NOT CHANGE WITHIN A GAME
READ (5,5) (GMNAMI I ) , 1=1,3)
5 FORMAT (3A8)
READ (5,6) LAND, MOLD, NONSTK, KAREA, I PLNT
6 FORMAT (2014)
IF (KAREA .EO. 0) GO TO 8
READ IN EQUAL AREAS IN ALL AGE CLASSES
NDX = MOLD + 1
DO 7 1=1, NDX
7 IACRE(I) = KAREA
ADJUST NUMBER OF ACRES IN OLDEST CLASS IF TOTAL AREA NOT MULTIPLE
OF KAREA
KDIFF = LAND - NDX » KAREA - NONSTK
I ACRE ( NDX ) = I ACRE I NOX ) + KDIFF
IACRE(l) = IACRE(l) + NONSTK
GO TO 10
READ IN UNEQUAL AREAS IN AGE CLASSES
8 READ (5,9) ( I ACRE ( I 1 ,1 = 1,180)
9 FORMAT (1514)
READ LIMITATIONS ON CUT
10 READ (5,11) (PRIDIVI I I
11 FORMAT (10F8.3)
18 -
READ (5,9) (MALCUTI I ) , 1=1, 10)
READ (5,11) ( FMRCHOI I ) , 1=1, 10)
READ (5,11) SHE LT,R ATE, CPLT,CTHN,CLOSS, ACCST, CUCST, BFCST, GROW, F I NL
READ (5,12) DEFOR, ANUL
12 FORMAT (2F8.5)
DEFGR1 = DEFOR * 100.0
PRINT CONDITIONS OF SIMULATIONS
WRITE (6,13)
13 FORMAT ( 1H1 ,//,46X,26HALTERNATIVES FOR THIS GAME)
WRITE (6,14) ( BATCH! I ), 1=1, 3)
14 FORMAT (1H ,45X,7HRATCH ,3A8)
WRITF (6,15) ITEST
15 FORMAT (1H , 45X , 4HTE S T , I 4 )
WRITE (6,16) ( GMNAM (I), 1=1, 3)
16 FORMAT ( 1H ,45X,6HGAME ,3A8)
WRITE (6,17) (DESCRl I ) , 1 = 1 ,5)
17 FORMAT (1H , 4 5 X , 5 A 8 , / / / / )
WRITE (6,18) NOYRS
18 FORMAT (1H , 45X , 24HNUMBER OF YEARS PER GAME, 14,//// )
WRITF (6,19) ( PRIDIV( I ) , 1=1, 10)
19 FORMAT (1H .15HCRITICAL PR I C E S , 1 2 X , 1 0F9 . 2 )
WRITE (6,201 (MALCUTI II , 1=1, 10)
20 FORMAT ( 1H , 13HALL0WABLE C UT , 1 1 X , 10 I 9 )
WRITF (6,21) (FMRCHDI I) , 1=1, 10)
?1 FORMAT ( 1H , 19HM I N I MUM CUTTING AGE , 5X , 1 OF 9. 0 , / / / / )
WRITE (6,22) LAND
22 FORMAT (1H , 23HACRES IN WORKING CIRCLE, 13X, I4,25X,27HC0STS IN FIRS
IT YEAR OF GAME)
WRITE (6,23) ACCST
23 FORMAT (1H ,69X,17HPER ACRE (ANNUAL), 8X,F9. 2)
WRITE (6,24) CUCST
24 FORMAT (1H , 38HM I N I MUM VALUES FOR INCLUSION IN TOTAL S , 3 1 X , 2 5HPER 1
100 CU. FT. HARVESTED, F9. 2)
WRITE (6,25) AGMRCH.BFCST
25 FORMAT ( 1H ,4X,22HAGE, FOR GROWING STOCK , 1 1 X , F3 . 0 , 2 9X , 1 3HPER M BD.
1 FT. , 12X.F9.2)
WRITE (6,26) BFMRCH, CTHN
26 FORMAT (1H ,4X,28HM BD. FT., FOR GROWING STOCK , 5X , F5 . 1 , 27X , 1 3HTH I N
1 ONE ACRE,12X,F9.2)
WRITE (6,27) COMCU.CPLT
27 FORMAT ( 1H ,4X,27HCU. FT., FOR COMMERCIAL CUT , 5X , F4. 0 , 29X , 1 4HPL A NT
1 ONE ACRE,11X,F9.2)
WRITE (6,28) COMBF.CLOSS
28 FORMAT (1H ,4X,29HM BD. FT., FOR COMMERCIAL CUT , 4X , F 5 . 1 , 27X , 19HC LE
1ANUP OF ONE ACRE,6X,F9.2)
WRITE (6,29) BF S ALV , RATE
29 FURMAT ( 1H .4X.22HM BD. FT., FOR S AL V AGE , 1 1 X , F 5 . 1 , 23X , 2 5HR ATE OF I
INCREASE IN COS T S , 4X , F9 . 2 , / / )
WRITE (6,30) IPLNT
30 FORMAT (1H , 22HACRES PLANTED ANNUALL Y , 1 4X , 1 4 , 2 5X , 3 5HREL AT I VE VALUE
1 OF INTERMEDIATE CUTS)
WRITE (6,31) DEFORl.CFPCT
31 FORMAT (1H , 30HPERCENT OF ACRES LOST ANNU ALL Y , 6X , F8 . 3 , 2 5X , 2 3HSTUMP
1 AGE PRICE, CD. FT.,2X,F9.2)
WRITE (6,32) SHELT.BFPCT
32 FORMAT ( 1H ,24HM BD. FT. IN SHELT ERWOOD , 1 3X , F5 . 1 , 27X , 2 3HSTUMP AGE P
lRICh, BD. FT.,2X,F9.2)
WRITE (6,33) FINL
33 FORMAT (1H , 19HREGENERAT ION PER IOD, 18X , F 5 . 1 , // )
WRITE (6,34) ANUL
34 FORMAT (1H , 29HPSEUD0RAND0M NUMBER GENER ATOR , 5 X , F 8 . 1 )
WRITE (6,35) GNTR
35 FORMAT (1H ,34X,F8.1)
RETURN
END
SUBROUTINE AREAS
COMMON ACCST, AGEO, AGMRCH, AN BDF (181),ANCUV(181),ANNET,BATCH(3),BDFC
1( 180) ,BDFO( 180) ,BFCST,CFMC( 180) ,CFMO( 1 80 ) , C S T AC , C ST VL , CUC ST , CUT AGE
2, DBHO.OENO, DESCRl 5 ) , DLE V, FMRCHD ( 10) ,GMNAM( 3) , G S VA LB , GS V ALC , G VLBF ,
3GVLCU,IACRE(180),IALCUT,IGAME,ISUM(18), ITEST, IVAR(26, 150), IY EAR,
4K0L (6) .LAND, LAST, MALCUT( 1 0 ) , MOL D, NACRE ( 1 80 ) , NG AME , NK OL S , N 1 , NONS TK ,
5N0YRS,PRET,PRIBD( 150) ,PRICF(150),PRIDIV(10) , RE TRN , ROT A , R ATE , GROW ,
6SITE,SUMM(6,25, 10) , TCOS T , SHE LT , YRLOS , F I NL , CLOS S , CTHN , CPLT , CFPC T ,
7VAR( 14, 150) ,VLBF,VLCU, AC AGE ( 1 000 ) , C YCL , KOUNT , 0 EFOR , ANUL , MI X , I PLNT ,
8RINT, THIN, BFMRCH, BFSALV, CO MCU,C0MBF,8FPCT, GNTR
GSVALB = 0.0
GSVALC = 0.0
GVLBF = 0.0
GVLCU = 0.0
DO 1 1=1,18
1 ISUMt I ) = 0
DO 2 1=1,1000
2 ACAGE ( I ) = 0.0
CONVERT ACRES IN EACH IACRE(I) TO INDIVIDUAL ACRES
JK = 0
DO 30 J=l,180
IF (JK . GE . LAND) GO TO 4
IF II ACRE ( J ) .LE. 0) GO TO 30
KL = JK + 1
JK = JK + IACRE(J)
DO 3 I=KL,JK
NAC = LAND +1-1
ACAGE ( NAC ) = J - 1
3 CONTINUE
30 CONTINUE
GET DISTRIBUTION OF ACRES BY AGE.
CHECK THAT NO ACRE IS OLDER THAN 179 YEARS
4 DO 7 K=l , LAND
IF ( ACAGE ( K ) .LE. 179.0) GO TO 6
WRITE (6,5)
5 FURMAT ( 1H1 ,47X,38HY0U WENT BEYOND AGE LIMIT OF 179 YEARS)
I YEAR = NOYRS - 1
RETURN
LM = AC AGE ( K ) + 1.0
N ACRE ( L M I = NACRE ( LM ) + 1
COMPUTE TOTAL ACREAGE BY 10-YEAR AGE CLASSES
DO 8 1=1,18
DO 8 J = l, 10
NS = 10 » (I - 11 + J
ISUM(I) = ISUM(I) + NACRE ( NS )
DO 10 M= 1 , L AND
IF ( ACAGE ( M ) .LT. AGMRCH) GO TO 10
I AG = ACAGE ( M ) + 1.0
IF ( ANBDF ( I AG ) .GE. BFMRCH) GO TO 9
GVLCU = GVLCU + ANCUV(IAG)
GO TO 10
9 GVLBF = GVLBF + ANBDF ( I AG )
10 CONTINUE
COMPUTE INITIAL NON-ZERO VALUES FOR 0UTPUT2
IVAR17, 1 )
- IVAR(8,1)
= GVLCU + 0.
= GVLBF + 0.
= IVAR(5,1)
) = I VARI6, 1 ]
) = NONSTK
= PRICF(l)
= PRIBDI 1 )
GVLBF * (PRIBD(l) - BFCST)
(GVLCU/100.0) » (PRICF(l) -
= GSVALC + GSVALB
= VARI 13, I ) + VAR( 12, 1 )
IVARI7.1
IVARI8, 1
IVARI9, 1
I VAR ( 10 ,
IVARI 11,
VAR( 1 ,1 )
VAR(2,1)
GSVALB =
GSVALC =
VAR( 13, 1
VAR ( 14, 1
DO 11 1=1,14
N = I + 11
IVAR(N.l) = ISUM(I)
IVAR(26,1) = ISUMI15
WRITE HEADINGS FOR TABLE OF INITIAL DISTRIBUTION OF ACRES BY AGE
WRITE (6,12)
12 FORMAT ( 1H1,////,38X,36HINITIAL DISTRIBUTION OF ACRES BY AGE)
WRITE (6,13) (BATCH! I ), 1=1 ,3)
13 FORMAT ( 1H ,45X,7HBATCH ,3A8)
WRITE (6,14) ITEST
14 FORMAT (1H ,45X,4HTEST, 14)
WRITE (6,15) ( GMNAM ( I ) , I = 1 , 3 )
15 FORMAT ( 1H ,45X,6HGAME ,3A8)
WRITE (6,16) (DESCRl I ), 1 = 1,5)
16 FORMAT (1H ,45X,5A8)
WRITE (6,17)
17 FORMAT (1H ,45X,16HYEAR WITHIN G AME , 3X , 1 HO , / // )
WRITE (6,18)
18 FORMAT (1H , 55X , 9HAGE ( YE AR ) )
WRITE (6,19)
19 FORMAT (1H , 4X, 1 1HAGE ( DECADE ), 7X , 1H0, 7X , 1H1 , 7X , 1H2.7X, 1H3,7X,1H4, 7
IX, 1H5,7X,1H6,7X,1H7,7X, 1H8.7X, 1H9 , 1 OX , 5HT0T AL , //)
WRITE NUMBER OF ACRES IN EACH
10-YEAR CLASSES
-YEAR AGE CLASS AND TOTALS OF
DO 21 J=l
IK = J -
NN = 10 * IK + 1
WRITE (6,20) I K , NACRE ( NN ) , NACRE (NN+1) , NACRE (NN+2) , NACRE ( NN + 3 ) , NACR
1EINN+4) , NACRE (NN+5) , NACRE ( NN+ 6 ), NAC RE (NN+7), NACRE (NN+ 8), NACRE (NN+9
2) , I SUM ( J )
20 FORMAT ( 1H , I 1 1 , 5X , 10 I 8 , I 1 5 , / )
21 CONTINUE
RETURN
END
SUBROUTINE YEARS
COMMON ACCST, AGEO, AGMRCH, ANBDF ( 1 8 1 ) , ANCU V ( 1 8 1 ) , ANNET , BATCH ( 3 ) , BDFC
1(180) ,BDFO( 180) , BFCST, CFMC (180) ,CFMO( 180) , C ST AC , C ST VL , CUCST , CUT AGE
2,DBH0,DEN0,0ESCR(5) , DL EV , FMRCHD ( 1 0 ) , GMN AM ( 3 ) , GS V AL B , G S VALC , G VL BF ,
3GVLCU, I ACRE ( 180 ) , IALCUT, I GAME , I SUM ( 18) , I TE ST , I VAR ( 26 , 1 50 ) , I YEAR,
4K0L ( 6) , LAND , LAST, MALCUTI 10), MOLD, NACRE ( 180 ),NGAME,NK0LS,N1, NONSTK,
5 NOYRS, PRET, PRIBDI 1 50 ) , PR I CF ( 1 50 ) , PR I D I V I 10 ) , RE TRN , ROT A , RA TE , GROW ,
6SITE,SUMM(6,25, 10) , TCOS T , SHE LT , YRLOS , F I NL , CLOS S , C THN , C PLT , CFPCT ,
7VAR(14,150),VLBF,VLCU,ACAGE(1000) ,CYCL, KOUNT, DEFOR, ANUL, MIX, IPLNT,
8RINT, THIN, BFMRCH, BFSALV, COMCU,COMBF,BFPCT, GNTR
GVLBF = 0.0
GVLCU = 0.0
LOSS = 0
NPLNT = 0
RETHV = 0.0
RETTH = 0.0
SCLOSS = 0.0
SCPLT = 0.0
SCTHN = 0.0
VLBF = 0.0
VLCU = 0.0
JCYCL = CYCL
IYRM = I YEAR + 1
MAKE
SCHEDULED ANNUAL PLANTING
IF (NONSTK .EQ. 0) GO TO 1
NPLNT = IPLNT
IF (NPLNT .GT. NONSTK) NPLNT =
NONSTK = NONSTK —NPLNT
APLT = NPLNT
SCPLT = APLT * CPLT
IF (NONSTK .EQ. 0) GO TO 1
KIM = LAND - NONSTK + 1 + LAST
19
IF (KIM .GT. LAND) KIM = KIM - LAND
ICH = KIM + NONSTK - 1
C
C AGE ZERO OF AREAS SUBROUTINE GIVES CONDITIONS AT END OF YEAR ZERO,
C SO INCREASE AGES ONE YEAR TO START SIMULATION
C
1 DO 2 1=1, LAND
ACAGE ( I ) = AC AGE ( I ) + 1.0
2 CONTINUE
IF (NONSTK . EO . 0) GO TO 4
C
C SUPPRESS AGE INCREASE FOR NONSTOCKED ACRES
C
DO 3 [=KIH, ICH
ACAGE ( I ) = 0.0
3 CONTINUE
4 IF (DEFOR .EO. 0.0) GO TO 15
C
C DETERMINE AREA DEFORESTED ANNUALLY
C
AKDX = LAND - NONSTK
YRLOS = (AKDX * DEFOR) + YRLOS
IF (YRLOS .LT. 1.0) GO TO 15
C
C GENERATE PSEUDORANDOM NUMBER FOR AGE OF ACRE DESTROYED
C
5 NDIV = (17.0 * ANUL + 3.01/128.0
NULL = ANUL
NULL = (17 * NULL + 3) - 128 * NDIV
ANUL = NULL
C
C CHECK THAT AGE EXISTS AND IS BETWEEN ONE AND OLDEST CURRENT AGE
C
IF (ANUL . L E . 0.0) GO TO 5
IF (ANUL .GT. ACAGE ( KOUNT ) ) GO TO 5
DO 6 M= 1 , LAND
KACR = M
IF ( ACAGE ( M ) . EO. ANUL) GO TO 7
6 CONTINUE
GO TO 5
C
C SET LOSS TO REDUCE CURRENT ALL0WA8LE CUT
C
7 LOSS = LOSS + 1
NONSTK = NONSTK + 1
YRLOS = YRLOS - 1.0
C
C SALVAGE VOLUME IF NOT LESS THAN BFSALV
C NO SALVAGE OR CLEANUP IF AGE LESS THAN AGEO.
C
IF (NULL .LE. Nl) GO TO 9
IF (IYEAR . EO. 1) MTHN = FMRCHDll)
NULL = NULL + 1
KULL = NULL - 1
IF (KULL .LT. MTHN) GO TO 8
SALVB = ANBDF ( NULL ) + BDFC(KULL)
IF (SALVB .LT. BFSALV) SCLOSS = SCLOSS + CLOSS
IF (SALVB .LT. BFSALV) GO TO 9
VL1F = VLBF + SALVB
RETTH = RETTH + SALVB * IPRIBDl IYRM) » BFPCT)
GO TO 9
8 SALVB = ANBDF ( NULL )
IF (SALVB .LT. BFSALV) SCLOSS = SCLOSS + CLOSS
IF (SALVB .LT. BFSALV) GO TO 9
VLBF = VLBF + SALVB
RETTH = RETTH + SALVB » (PRIBD(IYRM) » BFPCT)
C
C RENUMBER ACRES TO PUT ACRE LOST AT END OF AGE SEQUENCE
C WITH AGE ZERO
C
9 IF I KACR .NE. KOUNT) GO TO 10
LAST = LAST + 1
KOUNT = KOUNT + 1
ACAGE(LAST) = 0.0
GO TO 15
10 LUB = LAST - 1
IF (KACR .LT. LAST) GO TO 13
MNO = LAND - KACR
DO 11 J=1,MN0
JSUB = KACR + J
ISUB = JSUB - 1
ACAGE(ISUB) = ACAGE(JSUB)
11 CONTINUE
AC AGE ( L AND ) = ACAGE ( 1 )
DO 12 K=1,LUB
KAN = K + 1
ACAGE ( K ) = ACAGE ( KAN )
12 CONTINUE
ACAGE ( L AST ) = 0.0
GO TO 15
13 DO 14 M = K ACR , LUB
MOL = M + 1
ACAGE ( M ) = ACAGE ( MOL )
14 CONTINUE
ACAGE ( LAST ) = 0.0
IF (YRLOS .GE. 1.0) GO TO 5
C
C PREPARE SUBTOTALS FOR CURRENT YEAR AND CHECK THAT NO ACRE
C IS OLDER THAN 179 YEARS
C
15 DO 16 K=l,180
NACRE (K) = 0
16 CONTINUE
DO 19 K= 1 , LAND
IF (ACAGE(K) .LE. 179.0) GO TO 18
WRITE (6,17)
17 FORMAT ( 1HI ,/////, 47X.38HY0U WENT BEYOND AGE LIMIT OF 179 YEARS)
IYEAR = NOYRS
GO TO 42
18 LM = ACAGE ( K ) + 1.0
NACRE ( LM ) = NACRE ( LM ) + 1
19 CONTINUE
C
C DETERMINE ALLOWABLE CUT BASED ON BD.FT. STUMPAGE PRICE
C
DO 20 J=l,10
NSUB = J
IF (PRIBD(IYRM) .LE. PRIDIV(J)) GO TO 21
20 CONTINUE
21 IALCUT = MALCUT ( NSUB ) - LOSS
CUTAGE = FMRCHD ( NSUB )
C
C COMPUTE THINNINGS FOR ANNUAL CUT
C
MX Y = 0
MAC = CUTAGE
DO 24 I=N1,MAC, JCYCL
VBTH = 0.0
VCTH = 0.0
IF (I .GE. MAC) GO TO 25
MR = I + 1
IF (BDFC(I) .LT. COMBF) GO TO 22
VLBF 1 = NACRE ( MR ) « BDFC(I)
VLBF = VLBF + VLBF1
VBTH = VLBF 1
RETTH = RETTH + VBTH * (PRIBD(IYRM) ♦ BFPCT)
MX Y = MX Y + 1
GO TO 24
22 IF (CFMC(I) .LT. COMCU) GO TO 23
VLCU1 = N ACRE ( MR ) * CFMC(I)
VLCU = VLCU + VLCU1
VCTH = VLCU1
RETTH = RETTH + VCTH/100.0 * (PRICF(IYRM) » CFPCT)
MX Y = MXY + 1
GO TO 24
23 MXY = MXY + 1
SCTHN = NACRE I MR ) * CTHN + SCTHN
24 CONTINUE
25 MTHN = Nl + MXY * JCYCL
IF I IALCUT .LE. 0) GO TO 33
C
C COMPUTE VOLUME OF ACRES HARVESTED
C
KYR = IYEAR + 1 + F I NL
SHWD = SHELT + (SHELT * GROW)
ISHWD = SHWD + 0.5
DO 32 1=1, IALCUT
VBHV = 0.0
VCHV = 0.0
IF (LAST .LT. LAND) GO TO 26
LAST = 0
26 LAST = LAST + 1
IF (ACAGE(LAST) .GE. CUTAGE) GO TO 27
LAST = LAST - 1
GO TO 33
27 M = AC AGE ( L AST )
K = M + 1
KOUNT = KOUNT + 1
I SAFE = LAND + 1
IF (KOUNT .GE. ISAFE) KOUNT = 1
IF ( M .LT. MTHN) GO TO 28
VLBF2 = ANBDF ( K ) + BDFC(M) - SHELT
IF (VLBF2 .LT. COMBF) GO TO 29
VLDF = VLBF + VLBF2
VBHV = VLBF2
RETHV = RETHV + VBHV * PR I BD ( IYRM)
IVAR(4,KYR) = IVAR(4,KYR) + ISHWD
VAR(3,KYR) = VARI3.KYR) + (SHWD * PR I BD ( KYR ) )
VA R ( 7 , K YR ) = VAR ( 7 , KYR ) + SHWD * (BFCST + BFCST « RATE * FIND
GO TO 31
28 VLBF2 = ANBDF ( K ) - SHELT
IF (VLBF2 .LT. COMBF) GO TO 29
VL p c VL^F + VLBF2
vb:.. - vLBi-^
ET -1 v = RETHV + VBHV ♦ PRIBDI IYRM)
VAR i 4 , KYR ) = IVAR(4,KYR) + ISHWD
;AR(3,KYR) = V AR ( 3 , KYR ) +■ (SHWD * PR I BD ( K YR ) )
VAR ( 7 , KYR ) = VAR ( 7 , KYR ) + SHWD » (BFCST + BFCST » RATE * FIND
GO TO 31
29 IF (M .LT. MTHN) GO TO 30
VLCU2 = ANCUV(K) + CFMC(M)
IF (VLCU2 .LT. COMCU) GO TO 32
VLCU = VLCU + VLCU2
VCHV = VLCU?
RETHV = RETHV + VCHV/100.0 * PRICF(IYRM)
GO TO 31
30 VLCU2 = ANCUV(K)
IF ( VLCU2 .LT. COMCU) GO TO 32
VLCU = VLCU + VLCU2
VCHV = VLCU2
RETHV = RETHV + VCHV/100.0 * PRICF(IYRM)
31 ACAGEI LAST) = 0.0
32 CONTINUE
C
C COMPUTE GROWING STOCK VOLUME. USE CU.FT. IF VOLUME IS LESS
C THAN BFMRCH MBF
C
33 DO 37 1=1, LAND
IF (AC.GE(I) .LT. AGMRCH) GO TO 37
I AG = ACAGEI I ) + 1.0
IBG = I AG - 1
IF (IBG .LT. MTHN) GO TO 34
GBL1 = ANBDF ( I AG I + BDFC(IBG)
IF (GBL1 .LT. BFMRCH) GO TO 35
GVLBF = GVLBF + GBL1
GO TO 37
- 20
34 GBL1 = ANBDFI I AG )
IF (GBL1 .LT. BFMRCH) GO TO 35
GVLBF = GVLBF + GBL1
GO TO 37
35 IF (IBG .LT. MTHN ) GO TO 36
GCL1 = ANCUV(IAG) + CFMC(IBG)
GVLCU = GVLCU + GCL1
GO TO 37
36 GCL1 = ANCUVI I AG )
GVLCU = GVLCU + GCL1
37 CONTINUE
PREPARE FOR NEW TOTALS AND SUBTOTALS
DO 38 K=l,180
38 NACRE ( K ) = 0
DO 39 1=1,18
39 ISUM(I) = 0
DO 40 K= 1 , L AND
L M = AC AGE ( K ) + 1.0
NACREILM) = NACREILM) + 1
40 CONTINUE
COMPUTE TOTAL ACREAGE BY 10-YEAR AGE CLASSES
DO 41 1=1,18
DO 41 J=l,10
REMAINING SUBROUTINES
ISUMI I ) =
1 CONTINUE
COMPUTE VOLUMES AN6 VALUES AT
TO ANUAL
END OF CURRENT YEAR FOR TRANSFER
42 RETRN = RETTH J
CSTAC = LAND *
CSTVL = CUCST <
TCOST = CSTAC H
ANNET = RETRN -
GSVALB = GVLBF
RETHV
ACCST + SCPLT + SCTHN +
(VLCU/100.0) + BFCST «
CSTVL
TCOST
♦ (PRIBDI IYRM) - BFCST)
SCLOSS
VLBF
GSVALC = (GVLCU/100.0) * (PRICF(IYRM)
INCREASE COSTS ANNUALLY, IF DESIRED
+ (ACCST » RATE)
+ (BFCST ♦ RATE)
+ (CLOSS ♦ RATE)
(CPLT * RATE)
(CTHN « RATE)
+ (CUCST » RATE)
ACCST = ACCST
BFCST = BFCST
CLOSS = CLOSS
CPLT = CPLT +
CTHN = CTHN +
CUCST = CUCST
RETURN
END
SUBROUTINE 0UTPUT1
COMMON ACCST, AGEO, AGMRCH, ANBDFI 1 8 1 ) , ANCU V ( 1 8 1 ) , ANNE T , B A TCH ( 3 ) , BDF C
1 ( 180) ,BDFO( 180) , BFCST, CFMCI 180) ,CFMO( 180) , CSTAC , C ST VL , CUCS T , CUT AGE
2,DRH0,DEN0,DESCR(5) ,DLEV,FMRCHD( 10) , GMN AM ( 3), GSVALB, GSVALC, GVLBF,
3GVLCU, I ACRE! 180 ) , I ALCUT , I GAME , ISUM( 18) , I TE ST , I VAR ( 26 , 1 50 ) , I YEAR,
4K0LI 6) , LAND, LAST, MALCUTI 10) .MOLD, NACRE ( 180) , NG 4ME , NKOL S , N 1 , NON STK ,
5N0YRS.PRET, PRIBDI 150),PRICF(150),PRIDIV(10) , RE TRN , ROT A , R A TE , GROW ,
6SITE,SUMM(6,25, 10) , TCOST , SHE LT , YRLOS , F I NL , C LOS S , CTHN , C PLT , CF PCT ,
7VAR(14,150) , VLBF, VLCU, AC AGE I 1000) , C YCL , KOUNT , DEFOR , ANUL , M I X , I PLNT ,
SRI NT, THIN, BFMRCH, BFSALV,COMCU,COMBF,BFPCT,GNTR
WRITE TABLE HEADINGS
WRITE (6,1)
1 FORMAT I 1H1,////,46X,28HDISTRIBUTI0N OF ACRES BY AGE)
WRITE (6,2) (BATCH (I) ,1=1,3)
2 FORMAT (1H ,45X,7HBATCH ,3A8)
WRITE (6,3) ITEST
3 FORMAT (1H , 45X , 4HT EST , I 4 )
WRITE (6,4) (GMNAMI I) , 1=1, 3)
4 FORMAT (1H ,45X,6HGAME ,3A8)
WRITE (6,5) (DESCRI I ), 1=1, 5)
5 FORMAT (1H ,45X,5A8)
WRITE (6,6) IYEAR
6 FORMAT (1H ,45X,16HYEAR WITHIN GAME, 16,//)
WRITE (6,7)
7 FORMAT ( 1H , 55X , 9HAGE ( YE AR ) )
WRITE (6,8)
8 FORMAT (1H ,4X,11HAGE(DECADE),7X,1H0,7X,1H1,7X,1H2,7X,1H3,7X,1H4,7
1X,1H5,7X,1H6,7X,1H7,7X,1H8,7X,1H9,10X,5HT0TAL,//)
WRITE NUMBER OF ACRES IN EACH 1-YEAR AGE CLASS AND TOTALS OF
10-YEAR CLASSES
DO 10 J=l,18
K = J - 1
NN = 10 * K + 1
WRITE (6,9) K , NACRE ( NN ) , NACRE ( NN+ 1 ), NACRE < NN+2 ), NACRE ( NN+3 ), NACRE <
1NN+4) , NACRE ( NN+5 ) , NACRE ( NN+6 ),NACRE(NN+7),NACRE(NN+8),NACRE(NN+9),
? I SUM ( J )
9 FORMAT (1H ,111, 5X, 1018, 115,/)
10 CONTINUE
RETURN
END
SUBROUTINE ANUAL
CUfMON ACCST, AGEO, AGMRCH, ANBDFI 18 1) , ANCUVI 1 8 1 ) , ANNET , BATCH ( 3 ) , 8DFC
1( 180) ,BDFO( 180) , BFCST, CFMCI 180) ,CFMO( 180 ), CSTAC , CSTVL , CUCST , CUT AGE
2,DBHn,DEN0,DESCR(5),DLEV,FMRCHD(10) , GMNAMI 3 ), GSVALB, GSVALC, GVLBF,
3GVLCU, I ACRE I 180) , I ALCUT , I GAME, I SUM I 18) , ITEST, IVARI 26, 150) , IYEAR,
4K0L (6) , LAND, LAST, MALCUTI 10) , MOLD, NACRE! 180) ,NGAME,NK0LS,N1 »NONSTK,
5N0YRS,PRET,PRIBD( 150) , PR I CF ( 150 ) , PR I D I V ( 10 ) , RETRN , ROT A , R ATE , GROW ,
6S I TE , SUMM (6, 25, 10), TCOST, SHELT, YRLOS, FINL, CLOSS, CTHN, CPLT, CFPCT,
7 VAR I 14,150) , VLBF, VLCU, AC AGE (1000) , C YCL , KOUNT , LIE FOR , ANUL , M I X , I PLNT,
HRI NT, THIN, BFMRCH, BFS AL V , COMCU , COMBF , BFPCT.GNTR
CONVERT VOLUME AND AREA VALUES TO SUBSCRIPTED VALUES FOR USE BY
K = IYEAR
J = IYEAR +
IVAR(l.J) =
IVARI2.J) =
IVARI3.J) =
IBFT = VLBF
IVARI4.J) =
IVAR(5,J) =
IVAR(6,J) =
IVARI7.J) =
IVAR(8,J) =
IVARI9.J) =
IVARIIO.J) =
IVARIll.J) =
DO 1 1=1,14
N = I + 11
I ALCUT
CUTAGE
VLCU + 0.5
+ 0.5
I VAR I 4 , J )
IVARI5.K)
IVARI6.K)
GVLCU + 0.
GVLBF + 0.
I VAR ( 5 , J )
IVARI6.J)
■ NONSTK
IBFT
IVARI 3,
IVARI4,
I VAR I 7 , J )
• IVARI8.J)
1 IVAR(N.J) = ISUMII)
IVARI26.J) = ISUMI15)
ENTER MONEY VALUES IN ARRAYS FOR REMAINING SUBROUTINES
VAR(1,J)
VAR I 2, J)
VAR(3,J)
VAR ( 4 , J )
VAR(5,J)
VAR I 6 , J )
VAR I 7 , J )
VAR I 8 , J )
VAR I 9 , J )
VARI 10, J )
VARIll,
VAKI 12,
VARI 13,
VARU4,
RETURN
END
SUBROUTINE 0UTPUT2
COMMON ACCST, AGEO, AGMRCH, ANBDFI 181), ANCUV (181), ANNET , BATCH I 3 ) , BDFC
11180) ,BDFO( 180) , BFCST, CFMCI 180 ) ,CFMO( 1 80 ) , C STAC , C STVL , CUC ST , CUT AGE
2 , DBHO ,DEN0,DESCR(5),DLEV,FMRCHD( 10) , GMNAMI 3) , GSVALB , GS VALC , GVLBF ,
3GVLCU, I ACRE I 180) , I ALCUT , I GAME , I SUM ( 18 ), I TEST , I VAR I 26 , 150), IYEAR,
4K0L (6 ) , LAND, LAST, MALCUTI 10), MOLD, NACRE (180) , NGAME , NKOL S , Nl , NONSTK ,
5N0YRS,PRET,PRIBD( 150) , PR I C F I 150 ) , PR I D I V I 1 0 ) , RETRN , ROT A , R ATE , GROW ,
6SI TE, SUMM I 6, 25, 10 ), TCOST, SHELT, YRLOS, FINL, CLOSS, CTHN, CPLT, CFPCT,
7VAR1 14, 150) ,VLBF, VLCU, AC AGE I 1000) , C YCL , KOUNT , DE F OR , ANUL , M I X , I PLNT ,
8R I NT, THIN, BFMRCH, BFS ALV, COMCU, COMBF, BFPCT.GNTR
PRINT FIRST PAGE
PRICFI J)
PRIBDI J )
VAR (3, J) ■
VARI4.K) •
CSTAC
V AR ( 6 , K ) •
VAR I 7 , J ) •
VARI8.K) -
TCOST
= VARI10.K
= ANNET
= VARI12.K
= GSVALC +
= VARI 12, J
VARI5, J)
CSTVL
VAR(7,J)
+ V AR ( 9 , J )
+ VAR ( 11, J)
GSVALB
+ VARI 13, J )
NOYRS + 1
1 DO 15
LINE =
2 IF I M
.LT. 40) GO TO 11
3A8)
WRITE (6,3) (BATCH! I ), 1=1,31
3 FORMAT ( 1H1,//,46X,7HBATCH
WRITE (6,4) ITEST
4 FORMAT ( 1H , 45X , 4HTE ST , I 4 )
WRITE (6,5) (GMNAMI I ), 1=1, 3)
5 FORMAT (1H ,45X,6HGAME ,3A8)
WRITE (6,6) (DESCRI I ), 1=1,5)
6 FORMAT (1H ,45X,5A8)
WRITE (6,7)
7 FORMAT (1H ,/)
WRITE (6,8)
8 FORMAT ( 1H , 1 2X , 9HALL0WA BL E , 5X , 7HCUTT I NG , 8X , 10H ACTUAL CUT,10X,9HCU
1MUL CUT,10X,9HGRSTK VOL , 1 2 X , 9HT0T AL VOL)
WRITE (6,9)
9 FORMAT (1H , 2X , 4H YE AR , 9X , 3HCUT , 10X , 3HAGE , 7X , 6HCU . FT . , 5 X , 3HMBF , 6X , 6
1HCU.FT.,5X,3HMBF,6X,6HCU.FT., 5X , 3HMBF , 6X , 6HCU. F T . , 5 X , 3HMBF )
WRITE (6,10)
10 FORMAT (1H , 1 5X , 3H ( 1 ) , 1 OX , 3H ( 2 ) , 8X , 3H ( 3 ) , 7X , 3H ( 4 ) , 7X, 3H ( 5 ) , 7X , 3H ( 6
1 ) ,7X, 3H( 7) ,7X,3H(8) ,7X,3H(9) ,6X,4H( 10) ,//)
11 WRITE (6,12) LINE, ( IVARI I, J) , 1 = 1,10)
12 FORMAT I 1H ,16,112,113,112,19,3(111,19))
IF I J .LE. 1 ) GO TO 13
M = M + 1
IF ILL .LT. 10) GO TO 14
13 WRITE (6,7)
LL = 0
14 LL = LL + 1
15 CONTINUE
PRINT SECOND PAGE
M = 40
DO 23 J=1,N
LINE = J - 1
IF (M .LT. 40) GO TO 19
M = 0
WRITE (6,3) (BATCH! I ), 1=1,3)
WRITE (6,4) I TEST
WRITE (6,5) (GMNAMI I ) ,1=1,3)
WRITE (6,6) (DESCRI I ), 1=1,5)
WRITE 16,7)
WRITE (6,16)
16 FORMAT (1H , 1 1 X , 3HN0N , 46X , 1 1H AGE CLASSES)
WRITE (6,17)
17 FORMAT (1H ,2X,4HYEAR,5X,3HSTK, 114H 0-9 10-19 20-29 30-39 40
1-49 50-59 60-69 70-79 80-89 90-99 100-109 110-119 120-129
2 130-139 140-179)
WRITE (6,18)
18 FORMAT (1H , 1 OX , 1 1 7H I 1 1 ) (12)- (13) (14) (15) (16) (17)
1 (18) (19) (20) (21) (22) (23) (24) (25)
21
2 (26),//)
19 WRITE (6,20) L I NE , ( I VAR ( I , J ) , I = 1 1 , 26 )
20 FORMAT (1H ,16,217,16,817,18,419)
IF (J .LE. 1) GO TO 21
M = M + 1
IF (LL .LT. 10) GO TO 22
21 WRITE (6,7)
LL = 0
22 LL = LL + 1
23 CONTINUE
PRINT THIRD PAGE
H = 40
00 31 J=1,N
LINE = J - 1
IF (M .LT. 40) GO TO 27
M = 0
WRITE (6,3) (BATCH! J ) ,1=1,3)
WRITE (6,4) ITEST
WRITE (6,5) (GMNAMI I ) , 1=1 ,3)
WRITE (6,6) (DESCR(I),I=1,5)
WRITE (6,7)
WRITE (6,24)
24 FORMAT (1H , 1RX , 14HSTUMPAGE PRICE, 9X, 15HSTUMPAGE I NCOME , 1 3X , 10HARE
1A COSTS, 15X,12H VOLUME COSTS)
WRITE (6,25)
25 FOkMAT (1H , 2X , 4HYE AR , 9X , 1 0H1 00 CU. F T ., 5X , 3HMBF , 6X , 6H ANNUAL , 5X , 9HC
1UMULATED.6X , 6H ANNU AL , 5 X , 9HCUMUL AT ED , 6X , 6HANNU AL , 5X , 9HCUMUL ATED )
WRITE (6,26)
26 FORMAT ( 1H , 1 8X , 4H ( 27 ) , 8 X , 4H I 28 ) ,6X,4H(29) ,9X,4H( 30) ,9X, 4H ( 31 ) , 9X ,
14H( 32) ,9X,4H(33 ) ,9X,4H(34) ,/ / )
27 WRITE (6,28) LINE , ( VAR ( I , J ) , I = 1 , 8 )
20 FORMAT (1H , 1 6 , F 1 6 . 2 , F 1 1 . 2 , F 1 2 . 0 , F 1 1 . 0 , 2 ( F 1 5 .0 , F 1 1 . 0 ) )
IF (J .LE. 1) GO TO 29
IF (LL .LT. 10) GO TO 30
29 WRITE (6,7)
PRINT FOURTH PAGE
M = 40
DO 39 J=1,N
LINE = J - 1
IF (M .LT. 40) GO TO 35
M = 0
WRITE (6,3) (BATCH(I),I=1,3)
WRITE (6,4) ITEST
WRITE (6,5) (GMNAM(I),I=1,3)
WRITE (6,6) (DESCR(I),I=1,5)
WRITE (6,7)
WRITE (6,32)
32 FORMAT ( 1H , 18X , 10HT0TAL COST , 17X , 10HNET I NCOME , 1 3X , 1 3HCURR ENT VAL
1UE,9X,5HT0TAL(
WRITE (6,33)
33 FORMAT ( 1H , 2 X , 4H YE AR , 8X , 6HANNU AL , 5X , 9HCUMUL AT ED , 7X , 6H ANNUAL , 5X , 9H
1CUMULATED,7X, 13HGR OWING STOCK, 7X,9HNET WORTH)
WRITE (6,34)
34 FORMAT ( 1H , 1 5X , 4H ( 35 ) , 9X , 4H ( 36 ) , 10X , 4H ( 37 ) , 9X , 4H ( 3 8 ) , 1 3X , 4H ( 39 ) , 1
14X,4H(40) ,//)
35 WRITE (6,36) LINE,(VAR( I, J), 1 = 9, 14)
36 FORMAT (1H , 1 6 , F 14 . 0 , F 1 2 . 0 , F 1 5 . 0 , F 1 2 . 0 , 2F 1 8 . 0 )
IF (J .LE. 1) GO TO 37
M = M + 1
IF (LL .LT. 10) GO TO 38
37 WRITE (6,7)
LL = 0
38 LL = LL + 1
39 CONTINUE
RETURN
END
SUBROUTINE WORTH
COMMON ACCST,AGEO,AGMRCH,ANBDF( 181) , ANCUVI 1 8 1 ) , ANNET , B ATCH ( 3 ) , BDFC
1 ( 180) ,BDFO( 180) ,BFCST,CFMC( 180 ) ,CFMO( 1 8 0 ) , C ST AC , C ST VL , C UC ST , CUT AGE
?,DBHO,DENO,DESCR( 5) , OLE V , F MRCHD ( 1 0 ) , GMN AM ( 3 ) , GS V AL B , GS V ALC , GVL BF ,
3GVLCU , I ACRE ( 180) , I AL CUT , I G AME , I SUM! 18) , I TE ST , I VAR ( 26 , 1 50 ) , I YEAR,
4K0LI6) , LAND, LAST, MALCUT( 10) .MOLD, NACRE ( 1 80 ) , NG AME , NKOL S , N 1 , NONSTK ,
5N0YRS,PRET,PRIBD( 150) ,PRICF( 150) ,PRI0IV( 10) ,RETRN, ROTA, RATE, GROW,
6SITE,SUMM(6,25, 1 0 ) , TCOS T , SHELT , YRLOS , F I NL , CLOSS , CTHN , CPLT , CF PCT ,
7VAR( 14, 150) , VLBF.VLCU, AC AGE ( 1000) , C YCL , KOUNT , DEFOR , ANUL , M I X , IPLMT,
8RINT,THIN,BFMRCH,BFSALV,COMCU,C0MBF,BFPCT,GNTR
DIMENSION CRATE (20), DISCI 20) ,01 SGI 20) ,01 SI (20) , PREV( 20) , PWTH(20),R
1ATI0I20, 150)
DO 1 1=1,20
• CRATE ( I ) = 0.0
DISCI I )
DISG( I )
DISI ( I )
PREVI I )
PWTHI I )
DO 2 t=:
= 0.0
= 0.0
= 0.0
= 0.0
= 0.0
,20
DO 2 J=l,150
RATIO! I , J) = 0.0
PRF PARE ARRAY OF ALTERNATIVE RATES
CRATE! 1 ) = 0.010
DO 3 1=1,19
3 CRATF(K) = CRATE! I ) + 0.005
PREPARE INTEREST TABLE FOR PERIOD NOYRS
DO 4 K=l, NOYRS
YRS = K
PROD = A LOG I FACTR ) • YRS
RATIOIJ.K) = EXPIPROD)
4 CONTINUE
DISCOUNT GROWING STOCK VALUE AT NOYRS
DO 5 L=l,20
KL = NOYRS + I
DISG(L) = VAR! 13, KLI/RATIOIL, NOYRS)
5 CONTINUE
DISCOUNT ANNUAL COSTS AND RETURNS
DO
1,20
N)
PRESC = 0.0
PRESI = 0.0
SPRSC = 0.0
SPRSI = 0.0
DO 6 N=l, NOYRS
I = N + 1
PRESC = VAR(9, I )/RATIO(M,N
PRE S I = VAR(3, I 1/RATIOI
SPRSC = SPRSC + PRESC
SPRSI = SPRSI + PRE S I
CONTINUE
DISKM) = SPRSI
DISC(M) = SPRSC
CONTINUE
OBTAIN PRESENT WORTH AT EACH RATE
DO 4 J=l,20
FACTR = 1.0
DO 8 IJ=1,20
PREVI IJ) = DISI ( IJ) + DISGI IJ )
PWTH(IJ) = PREVIIJ) - DISC(IJ) - VAR(13,1)
CRATE(IJ) = CRATE(IJ) • 100.0
8 CONTINUE
WRITE (6,9)
9 FORMAT I 1H1,////,53X,29HPRESENT WORTH AND RATE EARNED)
WRITE (6,10) (BATCH! I), 1 = 1, 3)
10 FORMAT (1H ,52X,7HBATCH , 3A8 )
WRITE (6,11) ITEST
11 FORMAT (1H ,52X,4HTEST, 14)
WRITE (6,12) (GMNAMI I ), 1 = 1,3)
12 FORMAT (1H ,52X,6HGAME ,3A8)
WRITE (6,13) (DESCRI I ), 1=1,5)
13 FORMAT ( 1H ,52X,5A8)
WRITE (6,14) NOYRS
14 FORMAT ( 1H , 52X , 1 5HYE AR S IN PERIOD, 15,//)
WRITE (6,15) VAR(13,1)
15 FORMAT (1H , 1 1 X , 33HV ALUE OF INITIAL GROWING STOCK — J , F 1 0. 2 , // )
WRITE (6,16)
16 FORMAT (1H , 57X , 38HVALUES DISCOUNTED TO PRESENT (DOLLARS),/)
WRITE (6,17)
17 FORMAT (1H , 1 1 X , 8HC0MP0UND , 1 5 X , 6HFUTURE , 34X , 5HST0CK , 36X , 3HNET )
WRITE (6,18)
18 FORMAT ! 1H , 13X , 4HR ATE , 1 6X , 7HGR0W I NG , 1 5 X , 3H ALL , 1 7 X , 4HPLUS , 1 6X , 3HAL
1L,15X,7HPRESENT)
WRITE (6,19)
19 FORMAT (1H , 1 1 X , 9H ( PERC ENT ), 1 4X , 5HST0CK , 14X , 7H I NCOMES , 1 3X , 7H INCOME
1S,14X,5HC0STS, 15X, 5HW0RTH, / )
DO 21 1=1,20
WRITE (6, 20) CRATE! I ) , D I SG ( I ) , D I S I ( I ) , PRE V ( I ) , D I SC ( I ) , P WTH ( I )
20 FORMAT (1H , 1 2X , F5 . 1 , 1 2X , 5 ( F 10 . 2 , 1 OX ) , / )
21 CONTINUE
RETURN
END
SUBROUTINE SUMRY
COMMON ACCST.AGEO, AGMRCH, ANBDFI 181 ) , ANCUVI 181) , ANNET , BATCH I 3 ), BDFC
II 180) ,BDF0(180) ,BFCST,CFMC( 180) ,CFMO( 180) , C ST AC , C ST VL , CUCST , CUT AGE
2,DBH0,DEN0,DESCR(5),DLEV,FMRCHD( 1 0 I , GMN AM ( 3 ) , G S VALB , GS V ALC , G VLBF ,
3GVLCU, IACRE! 180) , I ALCUT, I GAME, I SUM ( 18) , I TE ST , I V AR ( 26 , 1 50 ) , I YEAR,
4K0LI6) ,LAND,LAST,MALCUTI 10) , MOLD, NACRE! 1 80 ) , NG AME , NKOL S , N 1 , NONSTK ,
5N0YRS,PRET,PRIBD( 1 50 ) , PR I CF ( 1 50 ) , PR I D I V ( 10 ) , RE TRN , ROT A , RATE , GROW ,
6SITE,SUMM(6,25, 10) , TCOST , SHELT , YRLOS , F I N L , CLOS S , CTHN , CPLT , CFPCT ,
7VAR( 14, 150) ,VLBF, VLCU, AC AGE ( 1 000 ) , C YCL , KOUNT , DE FOR , ANUL , M I X , I PLNT ,
8RINT,THIN,BFMRCH,BFSALV,C0MCU,C0MBF,BFPCT,GNTR
CONVERT IVARII.J) AND VAR(I,J) TO SUMMfl.J.K) AT END OF EACH GAME
LIM = 10 + NOYRS/10
DO 4 I=1,NK0LS
DO 4 J=1,LIM
K = KOL I I )
IF (J .GT. 10) GO TO 1
JJ = J + 1
GO TO 2
1 JJ = 10 •
2 IF (K .GT
SUMMI I , J, IGAME)
GO TO 4
3 K = K - 26
SUMMI I , J, IGAME) = VAR(K,JJ)
4 CONTINUE
WRITE SUMMARY TABLES IF ALL GAMES FINISHED
IF (IGAME .LT. NG AME ) GO TO 17
WRITE PAGE HEADINGS. SEPARATE PAGE FOR EACH COLUMN OF 0UTPUT2
IDENTIFIED IN INPUT1
DO 16 I=1,NK0LS
WRITE (6,5)
5 FORMAT (1HI,////,46X,26HC0MPARIS0N OF ALTERNATIVES)
WRITE (6,6) ( BATCH! I ), 1=1, 3)
6 FORMAT (1H ,45X,7HBATCH , 3A8 )
J - 10) + 1
26) GO TO 3
IVARIK, JJ)
- 22 -
ITEST
,45X,4HTEST,
1DESCRI I ) , 1
,45X,5A8>
14)
13,///)
WRITE (6,7)
7 FORMAT (1H
WRITE (6,8)
8 FORMAT (1H ,
K = KOL(I)
WRITE (6,9) K
9 FORMAT ( 1H , 45X , 8HCOLUMN
WRITE (6,10)
10 FORMAT (1H , 5X , 4HYEAR , 6X , 6HGAME 1,6X,6HGAME 2,6X,6HGAME 3,6X,6HG
IE 4,6X,6HGAME 5,6X,6HGAME 6,6X,6HGAME 7 , 6X , 6HGAME 8,6X,6HGAME 9,
2.7MGAME 10,//)
WRITE SUMM(I,J,K) FOR EACH OF FIRST 10 YEARS AND FOR END OF
EACH DECADE
DO 16
IF I J
JJ = J
GO TO 12
11 J J = 10 » ( J - 10)
12 IF ( H . L T . 5) GO TO 14
WRITE (6,13)
13 FORMAT (1H ,//)
M = 0
14 WRITE (6,15) JJ, ( SUMM( I , J,L) ,
15 FORMAT (1H , I 9 , F 1 1 . 0 , 9F 1 2 . 0 )
16 M = M+l
17 RETURN
END
APPENDIX 2
Output of Test Problem
PAGE TYPE I
YIELDS PER ACRE OF MANAGED, EVEN-AGED STANDS OF PONDE ROSA PINE IN THE BLACK HILLS
SITE
INDEX 60,
20-YEAR
CUTTING CYCLE, DENSITY
LEVEL
100
ENTIRE
STAND BEFORE AND
AFTER THINNING
PERIODIC CUT
AND MORT AL I T>
STAND
AGE
(YEARS)
TREES
NO.
BASAL
AREA
SQ.FT.
AVERAGE
D.B.H.
IN.
AVERAGE
HEIGHT
FT.
TOTAL
VOLUME
CU.FT.
MERCHANT-
ABLE VOLUME
CU.FT.
SAWTIMBER
VOLUME
MBF
TREES
NO.
BASAL
AREA
SQ.FT.
TOTAL
VOLUME
CU.FT.
MERCHANT-
ABLE VOLUME
CU.FT.
30
30
1000
496
110
73
4.5
5.2
20
20
818
534
201
201
0.
0.
504
37
284
0
40
496
104
6.2
28
1146
684
0.
50
50
496
286
133
86
7.0
7.4
35
35
1866
1209
1381
954
0.
0.
210
47
657
427
60
286
108
8.3
41
1806
1556
1.790
70
70
286
200
129
99
9.1
9.5
47
47
2501
1930
2241
1756
3.880
3.510
86
30
571
485
80
200
116
10.3
52
2515
2341
5.990
90
90
200
141
132
100
11.0
11.4
57
57
3144
2407
2940
2255
8.990
7.560
59
32
737
685
100
141
115
12.2
60
2902
2730
10.590
110
no
141
102
130
100
13.0
13.4
63
63
3457
2701
3263
2555
13.650
11.020
39
30
756
708
120
102
114
14.3
66
3216
3055
13.960
130
130
102
76
127
100
15.1
15.5
69
69
3741
2998
3565
2863
16.980
13.910
26
27
743
702
140
76
112
16.4
71
3445
3297
16.740
150
76
123
17.2
73
3887
3723
19.510
SAWTIMBER
VOLUME
MBF
PAGE TYPE 2
GROWING STOCK OF MANAGED BLACK HILLS PONDEROS A PINE
SITE INDEX 60, 20-YEAR CUTTING CYCLE
DENSITY LEVEL- 120 AND 100
VOLUMES PRESENT PER ACRE AT ENO OF EACH YEAR
MERCHANTABLE CUBIC FEET
DECADE
0
1
2
3
4
5
6
7
8
9
0
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.
1
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.
2
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.
3
201.0
249.3
297.6
345.9
394.2
442.5
490.8
539.1
587.4
635.7
4
684.0
753.7
823.4
893.1
962.8
1032.5
1102.2
1171.9
1241.6
1311.3
5
954.0
1014.2
1074.4
1134.6
1194.8
1255.0
1315.2
1375.4
1435.6
1495.8
6
1556.0
1624.5
1693.0
1761.5
1830.0
1898.5
1967.0
2035.5
2104.0
2172.5
7
1756.0
1814.5
1873.0
1931.5
1990.0
2048.5
2107.0
2165.5
2224.0
2282.5
8
2341.0
2400.9
2460.8
2520.7
2580.6
2640.5
2700.4
2760.3
2820.2
2880.1
9
2255.0
2302.5
2350.0
2397.5
2445.0
2492.5
2540.0
2587.5
2635.0
2682.5
10
2730.0
2783.3
2836.6
2889.9
2943.2
2996.5
3049.8
3103.1
3156.4
3209.7
11
2555.0
2605.0
2655.0
2705.0
2755.0
2805.0
2855.0
2905.0
2955.0
3005.0
12
3055.0
3106.0
3157.0
3208.0
3259.0
3310.0
3361.0
3412.0
3463.0
3514.0
13
2863.0
2906.4
2949.8
2993.2
3036.6
3080.0
3123.4
3166.8
3210.2
3253.6
14
3297.0
3339.6
3382.2
3424.8
3467.4
3510.0
3552.6
3595.2
3637.8
3680.4
15
3723.0
- 23 -
PAGE TYPE 2 (continued)
THOUSANDS OP BOARD FEET
0.
0.
0.
0.
1.790
3.510
5.990
7.560
10.590
11.020
13.960
13.910
16.740
19.510
0.
0.
.179
1.999
3. 758
6.290
7.863
10.896
11.314
14.262
14.193
17.017
.358
2.208
4.006
6.590
8. 166
11.202
11.608
14.564
14.476
17.294
0.
0.
c.
.537
2.417
4.254
6.890
8.469
11.508
11.902
14.866
14.759
17.571
0.
0.
0.
.716
2.626
4.502
7.190
8. 772
11.814
12.196
15.168
15.042
17.848
0.
0.
0.
0.
0.
.895
2.835
4.750
7.490
9.075
12. 120
12.490
15.470
15.325
18.125
0.
0.
0.
0.
0.
1.074
3.044
4.998
7.790
9.378
12.426
12.784
15.772
15.608
18.402
1.253
3.253
5.246
8.090
9.681
12.732
13.078
16.074
15.891
18.679
0.
0.
0.
0.
0.
1.432
3.462
5. 494
8.390
9.984
13.038
13.372
16.376
16. 174
18.956
0.
0.
0.
1.611
3.671
5.742
8.690
10.287
13. 344
13.666
16.678
16.457
19.233
PAGE TYPE 3 alternatives for this game
BATCH TEST PROBLEM
TEST 1
GAME EQUAL AREAS CUT ANNUALLY
MANAGED i THINNED AGE 30
NUMBER OF YEARS PER GAME 30
CRITICAL PRICES 99.00 0. 0. 0. 0. 0. 0. 0. 0. 0.
ALLOWABLE CUT 7000000000
MINIMUM CUTTING AGE 130 000000000
ACRES IN WORKING CIRCLE 915
MINIMUM VALUES FOR INCLUSION IN TOTALS
AGE, FOR GROWING STOCK 40
M BD. FT. , FOR GROWING STOCK 1.5
CU. FT. , FOR COMMERCIAL CUT 400
M BD. FT., FOR COMMERCIAL CUT 3.0
M BD. FT., FOR SALVAGE 1.5
COSTS IN FIRST YEAR OF GAME
PER ACRE (ANNUAL) .20
PER 100 CU. FT. HARVESTED .05
PER M BD. FT. 1.56
THIN ONE ACRE 25.00
PLANT ONE ACRE 30.00
CLEANUP OF ONE ACRE 25.00
RATE OF INCREASE IN COSTS .01
ACRES PLANTED ANNUALLY
PERCENT OF ACRES LOST ANNUALLY
M BD. FT. IN SHELTERWOOD
REGENERATION PERIOD
1
.040
4.0
10.0
RELATIVE VALUE OF INTERMEDIATE CUTS
STUMPAGE PRICE, CU. FT. 1.00
STUMPAGE PRICE, BD. FT. .85
PSEUDORANDOM NUMBER GENERATOR
PAGE TYPE 4
INITIAL DISTRIBUTION OF ACRES BY AGE
BATCH TEST PROBLEM
TEST 1
GAME EOUAL AREAS CUT ANNUALLY
MANAGED, THINNED AGE 30
YEAR WITHIN GAME 0
AGE ( DECADE )
- 24 -
PAGE TYPE 4 (continued)
DISTRIBUTION OF ACRES BY AGE
BATCH TEST PROBLEM
TEST 1
GAME EQUAL AREAS CUT ANNUALLY
MANAGED, THINNED AGE 30
YEAR WITHIN GAME 30
AGE(DECADE)
PAGE TYPE 5 BATCH TEST PROBLEM
TEST I
GAME EQUAL AREAS CUT ANNUALLY
MANAGED, THINNED AGE 30
ALLOWABLE
CUT
( 1 )
ACTUAL CUT
CU.FT. MBF
(3) 14)
CUMUL CUT
CU.FT. MBF
(5) (6)
GRSTK VOL
CU.FT. MBF
(7) 18)
CU.FT.
19)
MBF
(10)
130
130
130
130
130
130
130
130
130
130
16135
16135
16135
16135
16135
16135
16135
16135
14719
16135
16135
32270
48405
64540
80675
96810
1 12945
129080
143799
159934
182
272
363
454
545
636
727
805
897
145108
145108
145108
145108
145108
145108
145108
145108
145108
145108
4439
4439
4444
4444
4443
4447
4447
4446
4468
4467
161243
177378
193513
209648
225783
241918
258053
274188
288907
305042
4530
4621
4716
4807
4897
4992
5083
5173
5273
5364
130
130
130
130
130
130
130
130
130
130
16135
16135
16135
16135
16135
16135
16135
16135
16135
15450
114
127
122
114
127
122
122
128
123
119
176069
192204
208339
224474
240609
256744
272879
289014
305149
320599
1011
1138
1260
1374
1501
1623
1745
1873
1996
2115
145108
145108
145108
145108
145108
145108
145108
145108
144424
144354
4483
4483
4482
4499
4498
4497
4506
4505
4504
4517
321 177
337312
353447
369582
385717
401852
417987
434122
449573
464953
5494
5621
5742
5873
5999
6120
6251
6378
6500
6632
130
130
130
130
130
130
130
130
130
130
16135
16135
15427
16135
16135
16135
16135
16135
15708
16135
123
128
128
123
117
128
123
118
129
124
336734
352869
368296
384431
400566
416701
432836
448971
464679
480814
2238
2366
2494
2617
2734
2862
2985
3103
3232
3356
144284
144215
144145
144075
144006
143936
143866
143797
144154
144094
4516
4519
4520
4518
4531
4530
4528
4540
4539
4537
481018
497084
512441
528506
544572
560637
576702
592768
608833
624908
6754
6885
7014
7135
7265
7392
7513
7643
7771
7893
- 25 -
PAGE TYPE 5 (continued)
BATCH TEST PROBLEM
TEST 1
GAME EQUAL AREAS CUT ANNUALLY
MANAGED, THINNED AGE 30
NON
STK 0-9 10-19
(11) (12) (13)
20-29
( 14)
50-59
(17)
AGE CLASSES
60-69 70-79 80-89
(18) (19) (20)
90-99
(21)
100-109
( 22)
110-119 120-129 130-139 140-179
(23) (24) (25) (26)
BATCH TEST PROBLEM
TEST 1
GAME EQUAL AREAS CUT ANNUALLY
MANAGED, THINNED AGE 30
STUMPAGE PRICE
100 CU.FT. MBF
(27) (281
STUMPAGE INCOME
ANNUAL CUMULATED
(29) (30)
AREA COSTS
CUMULATED
(32)
VOLUME COSTS
ANNUAL CUMULATED
(33) (34)
2.50
2.50
2.50
2.50
2.50
2.50
2.50
2.50
2.50
2.50
15.20
17.80
16.80
13.40
14. 10
17.40
11.80
11.10
12.20
12.90
1784
2021
1880
1625
1688
1952
1482
1418
1325
1586
1784
3805
5686
7310
8999
10950
12433
13851
15176
16762
388
392
396
400
404
408
412
384
361
424
388
780
1176
1575
1979
2387
2799
3183
3543
3968
150
151
151
155
156
158
160
162
140
165
150
301
452
607
763
921
1081
1242
1383
1548
2.50
2.50
2.50
2.50
2.50
2.50
2.50
2.50
2.50
2.50
8.30
9.00
10.90
13.90
13.10
11.90
12.70
15.70
13.60
1566
1469
1512
1661
2192
2021
1862
2041
2346
2016
18329
19798
21310
22971
25163
27184
29046
31087
33433
35449
395
433
403
407
446
416
420
460
428
433
4363
4796
5199
5607
6053
6468
6888
7348
7776
207
232
225
214
240
233
236
248
240
236
1755
1987
2213
2427
2666
2899
3135
3382
3622
3859
2.50
2.50
2.50
2.50
2.50
2.50
2.50
2.50
2.50
2.50
12.10
15.20
16. 10
16.70
19.60
18.50
14.70
15.50
17.10
13.00
1904
2339
2471
2480
2726
2805
2235
2253
2617
2027
37353
39692
42163
44642
47368
50173
52408
54661
57278
59305
473
441
483
419
455
498
464
468
513
511
8682
9123
9606
10025
10479
10977
11440
11909
12421
12933
246
258
261
254
245
270
262
255
278
270
4104
4363
4623
4877
5123
5392
5654
5909
6187
6458
- 26 -
PAGE TYPE 5 (continued) batch test problem
TEST 1
GAME EQUAL AREAS CUT ANNUALLY
MANAGED, THINNED AGE 30
TOTAL COST NET INCOME CURRENT VALUE TOTAL
YEAR ANNUAL CUMULATED ANNUAL CUMULATED GROWING STOCK NET WORTH
(35) (36) (37) (38) (39) (40)
1 538 538 1247 1247 64099 65346
2 543 1081 1477 2724 75570 78294
3 547 1628 1334 4058 71139 75197
4 555 2182 1070 5128 55955 61083
5 560 2742 1128 6256 58989 65246
6 565 3308 1386 7643 73641 81284
7 572 3880 910 8553 48657 57210
8 545 4425 873 9426 45464 54890
9 501 4926 824 10250 50507 60757
1C 590 5516 997 11247 53551 64798
11 540 6056 659 11905 41104 53009
12 602 6658 565 12470 32952 45422
13 574 7232 657 13127 36007 49134
14 557 7789 708 13835 44593 58428
15 621 8410 1065 14900 58001 72901
16 592 9002 1020 15920 54315 70235
17 589 9591 840 16760 48919 65679
18 640 10231 939 17699 52431 70130
19 611 10841 1246 18945 65832 84777
20 600 11442 920 19865 56440 76305
21 660 12102 867 20732 49565 70296
22 630 12731 1156 21887 63512 85399
23 673 13405 1212 23099 67506 90605
24 611 14016 1347 24446 70106 94552
25 628 14644 1384 25830 83350 109180
26 695 15339 1437 27267 78246 105513
27 663 16002 1114 28380 60917 89298
28 649 16651 1039 29420 64611 94031
29 716 17367 1279 30698 71762 102461
30 717 18084 904 31603 53038 84641
PAGE TYPE 6 PRESENT WORTH AND RATE EARNED
BATCH TEST PROBLEM
TEST 1
GAME EQUAL AREAS CUT ANNUALLY
MANAGED, THINNED AGE 30
YEARS IN PERIOD 30
VALUE OF INITIAL GROWING STOCK — $ 60992.01
VALUES DISCOUNTED TO PRESENT (DOLLARS)
COMPOUND FUTURE STOCK NET
RATE GROWING ALL PLUS ALL PRESENT
(PERCENT) STOCK INCOMES INCOMES COSTS WORTH
1.0 39350.04 50396.60 89746.64 15451.27 13303.36
1.5 33931.66 46618.55 80550.21 14330.44 5227.76
2.0 29280.69 43222.05 72502.74 13320.23 -1809.50
2.5 25285.44 40163.29 65448.73 12408.09 -7951.37
3.0 21850.92 37403.89 59254.81 11583.01 -13320.21
3.5 18896.26 34910.23 53806.49 10835.34 -18020.87
4.0 16352.58 32652.78 49005.36 10156.61 -22143.27
4.5 14161.12 30605.60 44766.72 9539.35 -25764.64
5.0 12271.78 28745.87 41017.65 8977.00 -28951.36
5.5 10641.74 27053.48 37695.22 8463.76 -31760.55
6.0 9234.44 25510.68 34745.12 7994.52 -34241.41
6.5 8018.59 24101.82 32120.41 7564.75 -36436.36
7.0 6967.44 22813.03 29780.47 7170.46 -38382.01
7.5 6058.04 21632.05 27690.10 6808.09 -40110.31
8.0 5270.77 20548.01 25818.77 6474.49 -41647.73
8.5 4588.75 19551.25 24140.00 6166.86 -43018.87
9.0 3997.53 18633.19 22630.72 5882.70 -44243.99
9.5 3484.68 17786.21 21270.89 5619.79 -45340.91
10.0 3039.53 17003.50 20043.02 5376.15 -46325.14
10.5 2652.88 16278.99 18931.87 5150.01 -47210.15
- 27 -
PAGE TYPE 3
ALTERNATIVES FOR THIS GAME
BATCH TEST PROBLEM
TEST 1
GAME VARY CUT WITH PRICE
MANAGED, THINNED AGE 30
NUMBER OF YEARS PER GAME 30
CRITICAL PRICES
ALLOWABLE CUT
MINIMUM CUTTING AGE
12.00
5
130
99.00
10
120
ACRES IN WORKING CIRCLE 915
MINIMUM VALUES FOR INCLUSION IN TOTALS
AGE, FOR GROWING STOCK 40
M BD. FT., FOR GROWING STOCK 1.5
CU. FT., FOR COMMERCIAL CUT 400
M BD. FT., FOR COMMERCIAL CUT 3.0
M BD. FT. , FOR SALVAGE 1.5
COSTS IN FIRST YEAR OF GAME
PER ACRE (ANNUAL) .20
PER 100 CU. FT. HARVESTED .05
PER M BD. FT. 1.56
THIN ONE ACRE 25.00
PLANT ONE ACRE 30.00
CLEANUP OF ONE ACRE 25.00
RATE OF INCREASE IN COSTS .01
ACRES PLANTED ANNUALLY
PERCENT OF ACRES LOST ANNUALLY
M BD. FT. IN SHELTERWOOD
REGENERATION PERIOD
RELATIVE VALUE OF INTERMEDIATE CUTS
STUMPAGE PRICE, CU. FT. 1.00
STUMPAGE PRICE, BD. FT. .85
PSEUDORANCOM NUMBER GENERATOR
21.0
2222.0
PAGE TYPE 4
INITIAL DISTRIBUTION OF ACRES BY AGE
BATCH TEST PROBLEM
TEST 1
GAME VARY CUT WITH PRICE
MANAGED, THINNED AGE 30
YEAR WITHIN GAME 0
AGE(DECADE)
AGE(YEAR)
V 5
- 28 -
PAGE TYPE 4 (continued)
GAME VARY CUT WITH PRICE
MANAGED, THINNED AGE 30
YEAR WITHIN GAME 30
AGE I DECADE )
AGE I YEAR )
V 5
PAGE TYPE 5 BATCH TEST PROBLEM
TEST 1
GAME VARY CUT WITH PRICE
MANAGED, THINNED AGE 30
ALL JWABLE
ACTUAL CUT
CUMUL CUT
CU.FT. MBF
(5) (6)
GRSTK VOL
CU.FT. >
(7) I
TOTAL VOL
CU.FT.
19)
MBF
(10)
120
120
120
130
130
120
130
130
130
130
16135
16135
16135
16135
16135
16135
16135
16135
14719
16135
129
128
126
0
78
129
16135
32270
48405
64540
80675
96810
112945
129080
143799
159934
129
257
383
383
461
590
655
720
798
890
145108
145108
145108
145108
145108
145108
145108
145108
145108
145108
4389
4339
4294
4410
4426
4380
4413
4446
4468
4467
161243
177378
193513
209648
225783
241918
258053
274188
288907
305042
4518
4596
4677
4793
4887
4970
5068
5166
5266
5357
130
130
130
130
130
130
130
130
120
130
16135
16135
16135
16135
16135
16135
16135
16135
16135
15450
103
116
112
54
125
140
87
120
166
121
176069
192204
208339
224474
240609
256744
272879
289014
305149
320599
993
1109
1221
1275
1400
1540
1627
1747
1913
2034
145108
145108
145108
145108
145108
145108
145108
145108
144424
144354
4517
4551
4585
4636
4636
4635
4679
4678
4625
4638
321177
337312
353447
369582
385717
401852
417987
434122
449573
464953
5510
5660
5806
5911
6036
6175
6306
6425
6538
6672
130
120
120
120
120
120
130
120
120
130
16135
16135
15427
16135
16135
16135
16135
16135
15708
16135
115
160
159
153
156
166
111
155
178
336734
352869
368296
384431
400566
416701
432836
448971
464679
480814
2149
2309
2468
2621
2777
2943
3054
3209
3387
3430
144284
144215
144145
144075
144006
143936
143866
143797
144154
144094
4637
4588
4537
4484
4446
4393
4392
4353
4302
4401
481018
497084
512441
528506
544572
560637
576702
592768
608833
624908
6786
6897
7005
7105
7223
7336
7446
7562
7689
7831
- 29 -
PAGE TYPE 5 (continued) batch test problem
TEST 1
GAME VARY CUT WITH PRICE
MANAGED, THINNED AGE 30
0-9
( 12 :
10-19
(13)
30-39
115)
50-59
(17)
AGE CLASSES
60-69 70-79 80-89
(18) (19) (20)
90-99
(21 )
140-179
( 26)
BATCH TEST PROBLEM
TEST 1
GAME VARY CUT WITH PRICE
MANAGED, THINNED AGE 30
STUMP AGE PRICE
100 CU.FT. MBF
(27) ( 28 )
STUMPAGE INCOME
ANNUAL CUMULATED
(29) (30)
AREA COSTS
ANNUAL CUMULATED
(31) (32)
VOLUME COSTS
ANNUAL CUMULATED
(33) (34)
2.50
2.50
2.50
2.50
2.50
2. 50
2.50
2.50
2.50
2.50
15.20
17.80
16.80
13.40
14.10
17.40
11.80
11.10
12.20
12.90
2363
2682
2489
403
1501
2614
1169
1124
1325
1586
2363
5044
7533
7936
9438
12052
13221
14345
15670
17256
392
396
400
404
408
412
384
361
424
388
780
1176
1575
1979
2387
2799
3183
3543
3968
209
210
209
135
220
116
117
140
165
209
419
628
636
771
991
1 107
1224
1364
1530
2.50
2.50
2.50
2.50
2.50
2.50
2.50
2. 50
2.50
2.50
10. 10
8. 30
9.00
10.90
13.90
13. 10
11.90
12. 70
15. 70
13.60
1462
1385
1424
991
2151
2255
1439
1945
3024
2039
18718
20103
21527
22518
24669
26924
28363
30308
33333
35372
395
433
403
407
446
416
420
460
428
433
4363
4796
5199
5607
6053
6468
189
215
208
105
235
265
171
234
321
239
1719
1933
2141
2246
2481
2746
2916
3150
3471
3710
2.50
2.50
2.50
2.50
2.5C
2.50
2.50
2.50
12.10
15.20
16.10
16.70
19.60
18.50
14.70
15.50
17.10
13.00
1802
2815
2962
2970
3484
3499
2049
2821
3470
978
37174
39989
42952
45922
49406
52905
54954
57774
61245
62223
473
441
483
419
455
498
464
468
513
511
9123
9606
10025
10479
10977
11440
11909
12421
12933
230
319
320
311
322
345
236
330
381
102
3940
4259
4579
4890
5212
5557
5794
6123
6504
6606
- 30 -
PAGE TYPE 5 (continued)
LAT'H TFSI PROBLEM
TEST 1
GAME VARY
MANAGED, THI
TOTAL COST
ANNUAL CUMULATE!"!
(35) 136)
MFT INCOME
>iN'JAL CUHUL ATEU
I3M 138)
CURRENT VALUE
CROWING "STOCK
(39)
TOTAL
NET WORTH
CO)
597
602
604
408
539
628
528
501
501
590
597
1 199
1803
2211
2750
3378
3905
4406
4908
5497
1765
2080
1885
-5
963
1986
641
623
824
997
1765
3845
5730
5725
6688
8674
9316
9938
10762
11759
63417
73946
68856
55558
58777
72580
48313
45464
50507
53551
65182
77792
74586
61283
65465
81255
57628
55403
61269
65310
495
557
530
512
625
596
543
645
691
603
5993
6550
7079
7592
8216
8813
9356
10001
10692
11296
1093
1046
587
970
1843
940
12200
12597
13070
13548
14641
15687
16273
17243
19086
20026
41389
33399
36751
45848
59666
55867
50660
54307
67500
57855
53589
45996
49821
59396
74307
71554
66934
71550
86586
77882
664
710
753
689
705
770
658
724
787
549
1 1960
12670
13423
14112
14817
15587
16245
16969
17756
18305
887
1710
1791
1934
2065
2056
1085
1533
1 794
23
20913
22624
24415
26348
28413
30469
31554
33087
34861
34904
50797
64429
67750
69601
81844
75995
59188
62090
68202
51552
71710
87052
92165
95949
1 10258
106464
90742
95177
103083
86456
PAGE TYPE 6
PRESENT WORTH AND RATE EARNED
BATCH TEST PROBLEM
TEST 1
GAME VARY CUT WITH PRICE
MANAGED, THINNED AGE 30
YEARS IN PERIOD 30
VALUE OF INITIAL GROWING STOCK-
VALUES DISCOUNTED TO PRESENT (DOLLARS)
COMPOUND
RATE
(PERCENT)
1.0
1.5
2.0
2.5
3.0
3.5
4.0
4.5
5.0
5.5
6.0
6.5
7.0
7.5
8.0
8.5
9.0
9.5
10.0
10.5
FUTURE
GROWING
STOCK
38247.39
32980.84
28460. 19
24576.89
21238.62
18366. 76
15894.35
13764.30
11927.90
10343.54
8975.67
7793.90
6772.20
5888.29
5123.07
4460. 16
3885.51
3387.03
2954.35
2578.54
ALL
INCOMES
52747.65
48739.60
45142.72
41909.46
38998. 17
36372.40
34000.11
31853.20
29906.94
28139.54
26531.81
25066.81
23729. 57
22506.84
21386.90
20359.34
19414.93
18545.45
17743.61
17002.89
STOCK
PLUS
INCOMES
90995.04
81720.44
73602.91
66486. 35
60236.79
54739. 15
49894.46
45617.51
41834.84
38483.08
35507.48
32860.71
30501.77
28395. 13
26509.97
24819.50
23300.44
21932.48
20697.96
19581.44
ALL
COSTS
15614.36
14470.52
13440.43
12511. 14
11671.31
10910.99
10221.45
9595.00
9024.86
8505.07
8030.34
7596. 04
7198.01
6832.62
6496. 61
6187.10
5901. 53
5637.60
5393.28
5166. 75
NET
PRESENT
WORTH
14388. 67
6257.91
-829.53
-7016.80
-12426.53
-17163.85
-21319.00
-24969.50
-28182.04
-31014.00
-33514.87
-35727. 34
-37688.26
-39429.50
-40978.65
-42359.61
-43593. 10
-44697. 13
-45687.33
-46577.33
- 31 -
PAGE TYPE 7 COMPARISON OF ALTERNATIVES
BATCH TEST PROBLEM
TEST 1
MANAGED, THINNED AGE 30
COLUMN 10
4530
4621
4716
4807
4897
4518
4596
4677
4793
4887
4992
5083
5173
5273
5364
4970
5068
5166
5266
5357
5364
6632
7893
5357
6672
7831
0
90
100
110
120
130
140
150
COMPARISON OF ALTERNATIVES
BATCH TEST PROBLEM
TEST 1
MANAGED, THINNED AGE 30
COLUMN 40
65346
78294
75197
61083
65246
65182
77792
74586
61283
65465
81284
57210
54890
60757
64798
81255
57628
55403
61269
65310
64798
76305
84641
65310
77882
86456
90
100
110
120
130
140
150
- 32 -
co T3
CD C
60 rd
d
;3 *
rQ h a)
<-> «$ rd rH
?2lS
<u co to "O t!
^ CD CO ™ s_
• « S 8 »
a .S
8 h id
< 33
m ->-i
u
2 5 n CO
u CO t>
60 60
g a -
8 § g
rr? H *"
£ §^
rH .rH CO
-t-> £
rd 3
X +3
Ch cn
a s
CO O
co »
8 t
Hi (!
60 rci .5
CO w
£0) *
rv 60
1 5 <S a
o o
h o
11 p «3
ro
c X c
<U P cd
+3 co
d d
o m
60 •
d Q
r
CO
h h c
P i co CO
rd " CD
„ cd CD CO d
™ r; T5 CO CO
£ 5 M 5 S
s
cd ^
" S
CD w
rd
fi CD
H
33 rd
0 rd
2? o
C CD
g> CO
" *
O CO
-M 0
d >
d *
CD W
60 +j
cd co
d <0
cd u
CD
> U
T3
c d
rd .5
- T3
co CO
60 N
- ^ ^
T3 UJ +j
CO "3 <0
• d
TJ rH
H P >-}
td w
ft O
Oh '33
rd _
^ d
60 5
H rd
cd 33
CO ^
u
- oo
CO vO
< "33
O C
« H
a) co
ft m
Cd P rt
rd rt -
» ^ r2
co co
CD CO ™
* s s
J h o
■g .s u
s I 3
w § .S
3 - o S
cn *h o
<u oi «
3Q ■ S*
ft °
5 60 rd g
vh d
> u g» S
§ - 1 s
- 1 2 fl
o <5 0 d
rH to -j 10
C rl O Jj
h cd
■< H o
co H n, co
rH rrj CO d
r5 CO ^ %
rd
rd ^
! 3
rS S
CO vD
CO rH
O TJ C-i co
J D O
J_